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ABSTRACT 

This  report  consists  of  three  parts.  The  first  is  the  analysis  for 
minimization  of  the  temperature  distortion  due  to  the  thermocouple  cavity. 
The  error  is  minimized  or  reduced  to  zero  by  optimizing  the  combination 
of  cavity  diameter  and  depth  and  the  thermocouple  transport  properties  and 
size.  The  second  is  the  refinement  of  the  presently  available  computer 
program  for  prediction  of  the  surface  temperature  and  heat  flux  at  the 
inner  surface  of  the  pipe  by  inverting  the  temperature  response  measured 
by  an  interior  probe  close  to  the  heated  surface.  The  refinement  is 
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achieved  by  using  the  double  precision  format  in  the  program  and  adapting 

' the  dimensionless  formulation.  The  third  is  to  study  the  inversion  solution 

for  a large  time  duration  of  a time  dependent  surface  heat  flux.  The 

. solution  is  obtained  by  the  method  of  Laplace  transformation  with  the 

V ' 

convolution  Integral. 
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FOREWORD 

I 

I 

This  is  the  final  report  of  the  project  entitled  "Heat  Transfer 
Analysis  for  Unsteady  High  Velocity  Pipe  Flow."  The  project  was  funded 
by  Army  Research  Office  Grant  DAA-G29-76-G-0123  for  a six  month  period 
from  January  1976  to  June  1976. 

! 

The  work  done  for  project  is  reported  here  in  three  parts.  The  first 
is  the  analysis  for  minimization  of  the  temperature  distortion  due  to  the 
thermocouple  cavity.  The  error  is  minimized  or  reduced  to  zero  by  opti- 
mizing a combination  of  cavity  diameter  and  depth  and  the  thermocouple 
material  and  size.  The  second  is  the  refinement  of  the  presently  available 
computer  programs  for  prediction  of  the  surface  temperature  and  heat  flux 
at  the  inner  surface  of  the  pipe  by  inverting  the  temperature  measured  by 
an  Interior  probe  close  the  the  heated  surface.  The  refinement  is  achieved 
by  using  the  double  precision  format  in  the  program  and  adapting  the  dimen- 
sionless formulation.  The  third  is  to  study  the  Inversion  solution  for  a 
large  time  duration  of  a time  dependent  surface  heat  flux.  The  solution  is 
obtained  by  the  method  of  Laplace  transform  and  the  convolution  Integral. 

Each  of  the  above  three  subjects  is  reported  as  a part  of  the  present  report. 

The  investigator  would  like  to  thank  Mr.  D.  M.  Thomsen  of  General 
Rodman  Laboratory,  Rock  Island  Arsenal,  Rock  Island,  Illinois  for  his  con- 
stant participation  in  the  course  of  the  research. 
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NOMENCLATURE 


specific  heat 
diameter  of  thermocouple 


diameter  of  cavity 


thickness  of  the  disk 


constant  heat  flux 


temperature 


real  time 


coordinate  along  and  normal  to  the  heated  surface 


Error  % [0  (at  cavity  base)  - 6 (at  edge  of  disk)]  x 100,  or  distorted 

temperature  divided  by  QD/tc 


Greek  letters 


distance  between  the  base  of  cavity  to  the  heated  surface 

dimensionless  temperature  Tk^^/QD 

dimensionless  time 

thermal  dlffusivity 

thermal  conductivity 

density 


subscripts 


1 the  disk  material 

2 Insulation  material 

3 thermocouple  material 

t thermocouple 
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ON  MINIMIZATION  OF  TPMPF.RATURF  DISTORTFON 
IK  THE  THERMOCOUPLE  CAVin 
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SUMMARY  - PART  I 


When  a thermocouple  Is  embedded  in  a solid  to  measure  temperature 
a distortion  of  temperature  is  created  because  of  the  thermocouple  cavity. 
For  a given  solid  under  measurement  the  temperature  distortion  at  the  base 
of  the  thermocouple  cavity  can  be  minimized  by  a proper  choice  of  the 
thermocouple  size  and  its  material,  the  cavity  diameter  and  the  depth  of 
the  cavity.  This  optimum  combination  is  solved  in  this  study  by  a finite 
element  analysis  for  the  case  that  the  thermocouple  cavity  is  drilled  into 
the  center  of  a disk.  The  disk,  initially  at  a uniform  temperature,  is 
then  heated  with  a constant  heat  flux  on  one  suface  and  is  insulated  on 
the  other  surface.  The  calculated  result  covers  a range  of  thermal  con- 
ductivity and  dlffusivlty  for  most  commonly  used  thermocouples.  For  ease 
in  practical  applications  a simple  formula  for  determining  the  optimum 
ratio  of  the  thermocouple  diameter  to  that  of  the  cavity  is  given  as  a 
function  of  density  - specific  heat  product  and  thermal  conductivity  of 
solid  material,  thermocouple,  and  insulating  material. 


LIST  OF  FIGURES  - PART  T 


PACE 


Figure  1 Geometric  Representation  of  Problem  19 

Figure  2 Finite  Element  Idealization  20 

Figure  3a  Thermal  Conductivity  of  Thermocouple  Materials  21 

Figure  3b  Variation  of  Thermal  Conductivity  Ratio  and  Density  - 22 

Specific  Heat  Ratio 

Figure  4 Temperature  Distribution  with  Insulation  Material  Filled  23 

the  Cavity 

Figure  5 Temperature  Distribution  with  Thermocouple  Material  Filled  24 

the  Cavity 

Figure  6 Temperature  Distribution  with  Thermocouple  Material  Partially  25 

Filled  the  Cavity 

Figure  7 Percentage  Error  vs.  d^/d  Ratio  for  Various  and  pc  26 

Ratios  e/D  = 0.1 

Figure  8 Percentage  Error  vs.  d^/d  Ratio  for  Various 
Ratios  e/D  = 0.04 

Figure  9 Percentage  Error  vs.  d(./d  Ratio  for  Various 
Ratios  e/D  = 0.02 


LIST  OF  TABLES 


PART  I 


PAGE 


Table  la  Commonly  Used  Thermocouples 

Table  lb  Thermal  Properties  of  Thermocouple  Materials 

Table  2 Percentage  Error  of  Temperature  at  Cavity  Base  c/D  = 0.1 

Table  3 Percentage  Error  of  Temperature  at  Cavity  Base  c/D  = O.OA 


Table  4 Percentage  Error  of  Temperature  at  Cavity  Base  e/D  - 0.02 


29 

30 

31 
34 
37 


INTRODUCTION 


A direct  measurement  of  transient  surface  temperature  and  heat  flux 
Is  often  difficult.  For  example,  a surface  Involves  two  modes  of  heat 
transfer,  say,  radiative  and  convective  heat  transfer.  In  this  case  If 
the  measuring  probe  has  a different  radiative  property  than  that  of  the 
surface,  erroneous  measurements  will  result.  A piston  or  projectile 
sliding  over  the  cylindrical  surface  Is  another  case  where  the  direct 
measurement  at  the  surface  Is  difficult.  A surface  Involving  melting  or 
ablation  Is  also  difficult  to  make  direct  measurement.  Therefore,  In- 
direct estimation  by  Inverting  the  temperature  history  Inside  the  heat 
conducting  solid  as  measured  by  a thermocouple  Is  often  used  for  pre- 
diction of  the  surface  temperature  and  heat  flux.  Beck  [1],  Hernnlng 
and  Parker  [2],  Frank  [3],  Imber  and  Khan  [4],  Stolz  [5],  Chen  and 
Thomsen  [6]  have  developed  different  Inversion  solutions  for  this  pur- 
pose. All  of  these  solutions  assumed  that  the  cavity  drilled  into  the 
solid  does  not  distort  the  true  temperature  distribution.  Thus,  it  Is 
important  that  the  temperature  measurement  by  an  Interior  probe  is 
accurate  and  involves  least  distortion  or  error.  Theoretically  Beck  [7], 
Masters  and  Stein  [8],  Burnett  [9],  and  Chen  and  Li  [10]  studied  the  distor- 
tion of  the  temperature  field  In  the  presence  of  a thermocouple  and  its 
cavity.  Experimentally  Chen  and  Danh  [11]  showed  that  appreciable  dis- 
tortion, say  10%,  of  temperature  field  may  exist  for  a normal  Implant 
of  the  thermocouple  Into  a solid  body.  From  studies  of  Chen  and  LI  [10] 
and  Beck  [7]  they  found  that  with  a proper  combination  of  the  thermocouple 
cavity  diameter.  Its  depth,  and  the  thermocouple  material  and  its  diameter. 


the  distortion  of  temperature  field  with  respect  to  space  or  time  can 
he  minimized  if  not  eliminated.  In  this  report  we  study  the  optimum 
combination  of  these  parameters  such  that  at  a given  situation  one 
knows  what  is  the  best  combination  to  use  and  what  is  the  magnitude  of 
the  temperature  distortion. 


FORMULATION  OF  PROBLEM 

In  the  present  study  we  consider  a disk  depicted  in  Figure  1 which 
has  a thickness  D and  is  drilled  a cavity  of  a diameter  to  a depth 
of  c distance  from  the  heating  surface.  The  heat  flux  Q is  assumed  to 
be  constant  and  the  upper  surface  of  the  disk  is  assumed  to  be  insulated. 

A thermocouple  of  a diameter  d^  is  then  welded  on  the  cavity  base. 
Furthermore,  the  disk  may  be  thought  to  approximate  a cut  out  from  a 
hollow  cylinder  if  the  radius  of  the  disk  is  small  compared  with  the  radius 
of  the  cylinder.  The  diameter  of  the  disk  is  choosen  to  be  D where  the 
temperature  distortion  due  to  the  thermocouple  cavity  becomes  negligible. 
For  this  to  be  true  one  needs  to  restrict  the  ratio  of  cavity  diameter 
to  the  disk  diameter  d/2D  be  small.  The  unfilled  cavity  can  be  air  or 
insulating  material. 

The  basic  idea  to  minimize  or  to  eliminate  the  temperature  distortion 
is  based  on  a proper  choice  of  the  thermocouple  size  and  the  material 
which  has  a higher  thermal  conductivity  than  that  of  the  disk  so  as  to 
conduct  more  heat  away  at  the  cavity  base  balancing  the  insulation  effect 
of  the  Insulator  in  the  cavity. 

Let  X and  Y be  respectively  the  coordinate  along  and  normal  to  the 
heated  disk  surface  and  the  Y a>is  coincide  with  the  axis  of  the  cavity. 
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for  the  thermocouple  (subscript  "3") 

39  “3  1 39  ^ 3^g\ 

2 

where  "t  = a^^t/D  Is  the  dimensionless  time,  x = X/D  the  dimensionless 
radial  coordinate,  y = Y/D,  the  dimensionless  distance  normal  to  the 
heated  surface,  a^,  and  are  respectively  the  thermal  dlffuslvity 

of  the  disk,  insulating  material,  and  the  thermocouple.  The  dimension- 
less temperature  9 is  defined  as  T<j^/QD  where  T is  the  temperature  above 
the  initial,  uniform  temperature. 

The  intlal  temperature  of  the  disk  is  then 

9 (x,  y,  0)  = 0 (A) 


a constant  heat  flux  at  the  lower  surface. 


y = 0 — 

ay 


the  Insulation  at  the  upper  surface. 


n 90 

y “ ° W 


the  zero  temperature  distortion  at  the  edge  of  the  disk. 


I 96 

X = 1 — 

3x 


and  the  axisymmetric  condition  at  the  cavity  axis 


n 90 
^ = 0 ^ 


The  condition  (7)  of  the  zero  temperature  distortion  was  verified  by 
Chen  and  LI  [10]  in  their  early  calculation  when  the  cavity  diameter  is 
one  tenth  of  the  disk  diameter.  In  addition  to  the  above  boundary 
condition  the  temperature  and  heat  flux  at  the  Interface  of  the  disk, 
thermocouple,  and  insulating  material  are  taken  to  be  continuous. 


ANALYSIS 


There  are  five  parameters  that  can  be  varied  for  the  present  anal- 
ysis. They  are  (a)  the  dimensionless  distance  from  the  base  of  the 
cavity  to  the  heated  surface  e/D,  (b)  the  size  of  the  cavity  d/D, 

(c)  the  ratio  of  the  thermocouple  diameter  to  that  of  the  cavity,  dj./d. 


! 
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(d)  the  thermal  conductivity  ratio  comes  from  the 

continuity  of  heat  flux  at  Interfaces.  (e)  the  r.itlo  of  the  product  of 
density  and  specific  heat  equivalently  to  the  ratio  of 

of  thermal  dlfflslvlty  “3/012  subscripts  1,  2 and  3 denote 

the  disk.  Insulation  and  thermocouple  materials. 

Because  of  the  complexity  of  the  geometry  and  the  multiplicity  of 
material  the  method  of  finite  element  technique  as  discussed  by  Wilson 
and  Nikel  [12]  is  adapted  with  the  aid  of  a computer  program  developed 
by  Wilson  [13].  The  present  problem  is  subdivided  into  finite  element 
as  required  by  the  method,  see  Figure  2.  The  dimensionless  pie  section 
Is  subdivided  into  121  finite  elements  with  12  dividing  lines  on  both 
coordinates.  Each  element  Is  defined  by  four  nodal  points  where  nodal 
points  are  denoted  by  intersections  of  the  dividing  lines  and  numbered 
as  shown  in  Figure  2.  The  material  property  corresponding  to  each 
element  is  then  assigned  to  the  program  developed  by  Wilson  [13]  . 

The  solution  at  each  nodal  with  respect  to  time  is  then  obtained. 

For  numerical  calculation  three  typical  values  of  the  distance 
from  heated  surface  to  the  base  of  the  cavity  e/D  are  chosen  to  be  0.04, 
0.1  and  0.2.  The  cavity  diameter  is  fixed  at  one  tenth  of  the  disk 
diameter.  The  thermocouple  to  cavity  diameter  ratio  d^/d  is  made  to 
vary  0,  0.2,  0.4,  0.6,  0.8  and  1.0.  Regarding  the  range  of  the  ratio 
of  the  thermal  conductivity  and  the  ratio  of  the  product  of  density  and 
specific  heat  we  surveyed  these  ratios  for  the  commonly  used  thermocouples 
In  Table  1 [14]  and  plotted  In  Figure  3.  Figure  3b  shows  the  ratio  of 
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thermocouple  coiuhicL J vlty  to  tluit  of  the  cUrR  versus  the  ratio  of 

density  “ sp.'clfie  heat  imnluct  f wliere  tlu’  value  of  the  conduc- 

tivity Is  taken  to  he  tlie  averaf5e  value  between  200  and  H00°K.  One  secs 
that  for  most  practical  situations  tlie  ratio  of  density  - specific  licat 
product  is  approximately  constant  at  0.7  except  when  the  con- 
ductivity ratio  is  small.  Thus  the  value  of  P3*^3/P3^]^ 

for  calculation  are  chosen,  as  shown  in  triangular  symbols  of  Figure  3b, 
to  cover  the  practical  range.  The  corresponding  value  for  and 

k^/k^  for  the  Insulation  material  are  chosen  to  be  fixed  at  0.5  and  0.005 
which  is  a typical  value  for  Teflon  Insulating  material  and  is  also  ap- 
proximately the  order  of  magnitude  for  air. 

RESULTS  AND  DISCUSSIONS 

Numerical  results  of  the  calculations  are  presented  in  Tables  2 to 
4 and  Figures  4 to  9.  The  percentage  error  of  temperature  is  defined 
as  the  distorted  temperature  divided  by  a reference  temperature  defined 
by  QD/Kj^.  Tables  2 to  4 give  the  percentage  error  of  temperature  dis- 
tortion) as  ;i  function  of  time  for  different  values  of  the  parameters 
d^/d  (0  to  1.0),  (0.5  to  10),  03^3/93^^2^ 

(0.02  to  0.1). 

Figures  4,  5 and  6 show  the  tlirec  typical  temperature  distributions 
in  the  steel  difik  near  the  thermocouple  junction  for  the  case  c/D  - 0.06 
d/2D  - 0.1  at  the  time  T = 0.08.  Figure  4 is  the  temperature  distribu- 
tion when  the  cavity  is  filled  entirely  with  the  insulation  material 
(Ttrflon  "2"  \..^h  ^ = 0.005,  2^  ^ 1 ~ ^ • '"0  • Tlie  dimensionless  isotherms 


10 


Tic  /QD  shows  the  distortion  of  temperature  distribution.  One  sees  that 
1 

the  Insulation  effect  on  the  heat  transfer  near  the  cavity  base  not  only 
creates  a much  higher  Junction  temperature  of  Tk^/QD  = 0.3A2  than  the 
undlstorted  one  of  0.255  at  the  edge  of  the  disk  giving  an  error  of  8.7% 
but  also  causes  a hot  spot  at  the  heating  surface  with  a higher  tempera- 
ture of  Tk^/QD  =■  0.374  over  the  undistorted  one  of  0.311.  On  the  other 
hand  the  temperature  distribution  in  the  Insulation  material  is  much 
lower  than  the  true  temperature  creating  a large  temperature  gradient 
at  the  base  of  the  cavity.  Figure  5,  contrary  to  Figure  4,  is  the  tem- 
perature distribution  when  the  cavity  is  completely  filled  with  thermo- 
couple material  whose  thermal  conductivity  is  ten  times  larger  than  the 
disk  material  e.g.,  copper  versus  steel.  Now  over-conduction  of  heat 
by  the  thermocouple  has  created  a cold  spot  at  the  base  of  the  cavity 
giving,  Tic^/QD  of  0.17  versus  the  undlstorted  one  of  0.255  with  an 
error  of  8.5%  as  well  as  at  the  heating  surface  with  Tic^/QD  of  0.236 
versus  0.311.  The  temperature  distribution  in  the  thermocouple  now 
becomes  higher  than  the  undlstorted  one  in  the  disk.  By  properly 
choosing  the  ratio  of  thermocouple  diameter  to  that  of  the  cavity  one 
may  minimize  these  distortions  of  the  temperature  response  at  the  base 
of  the  cavity.  This  is  shown  in  Figure  6 where  the  thermocouple 
diameter  d^  is  chosen  to  be  0.4  of  the  cavity  diameter  d with  *“  10 » 

e.g.  copper-steel  combination.  Figure  6 shows  that  distortions  at  the 
base  of  the  cavity  and  at  the  heating  surface  are  almost  eliminated  giving 
the  temperature  Tk^/QD  of  0.245  and  0.320  with  respectively  the  error  of 
1%  and  0.9%. 
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In  order  to  examine  the  details  of  the  distorted  temperature  response 
at  the  base  of  the  cavity  we  tabulated  the  results  in  Tables  2,  3 and  4 
and  plotted  the  error  percentage  as  function  of  the  ratio  of  the  thermo- 
couple diameter  to  that  of  the  cavity  for  different  cavity  depth,  time, 
and  thermal  conductivity  In  Figures  7,  8 and  9.  In  these  Figures  the  pc 
ratio  ranges  from  0.7  to  1.3  which  covers  most  of  the  practical  applica- 
tions . 

In  the  case  when  the  pc  ratio  is  equal  to  or  less  than  one  the  errors 
In  temperature  response  at  the  base  of  the  cavity  In  these  figures  are 
all  positive  or  overheat  for  ^ This  is  because  the  thermal 

conductivity  of  the  thermocouple  is  less  than  that  of  the  disk  and  the 
heat  capacity  pc  of  the  thermocouple  Is  also  small.  Therefore,  no  extra 
conduction  of  heat  can  be  achieved  by  the  thermocouple  to  compensate  for 
the  blocking  of  the  heat  transfer  by  the  insulation  material  In  the 
cavity.  On  the  other  hand,  if  ^ ^ error  of  temperature  varies 

from  positive  value  for  d^/d  ■ 0 to  some  negative  value  as  d^/d  approaches 
1.  Thus  for  > 1 a properly  chosen  combination  of  thermocouple 

and  Insulation  material  can  minimize  the  error.  For  example.  In  Figure  7 
a combination  of  ^ ~ d^/d  * 0.5 

produces  almost  negligible  error.  This  combination  shown  to  be  optimum 
at  T ■ 1.0  in  Figure  7 Is  also  optimum  for  other  time  periods  (see  Table 
2).  Therefore  once  an  optimum  combination  of  parameters  is  chosen  it  is 
valid  throughout  the  entire  transient  period  of  an  experiment.  From 
Figures  7,  8 and  9 one  can  also  see  that  the  optimum  ratio  of  d^/d  which 
gives  zero  temperature  error  decreases  as  the  increases. 
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This  implies  that  fpr  the  thermocouple  with  a larger  thermal  conductivity 
a smaller  diameter  Is  sufficient  to  eliminate  the  temperature  distortion. 


The  result  shown  in  Figures  7,  8 and  9 can  in  general  be  adapted  for  use 
In  practical  application  to  choose  the  size  of  thermocouple  and  cavity, 
the  thermocouple  material  and  the  depth  of  the  cavity  to  be  drilled. 

As  mentioned  earlier  that  when  ^ ^ ^ errors 


of  the  temperature  response  at  the  base  of  the  cavity  are  all  positive. 
However  we  found  (see  Tables  2.2,  3.2  and  4.2)  that  if  ratio 

is  made  large  enough  during  the  transient  period  the  error  of  temperature 
response  at  the  base  of  the  cavity  may  indeed  become  negative  even  when 
<3/^1  < 1.  Physically  although  the  thermal  conductivity  of  the  thermo- 
couple Kj  is  smaller  than  that  of  the  disk  material,  but  with  a larger 
heat  caplcltance  P3p3^P3‘^j^  ^ ^ thermocouple  is  still  capable  of 
absorbing  extra  heat  flux  and  hence  eliminates  the  temperature  distor- 
tion at  the  cavity  base  during  the  transient  period.  To  illustrate  this 
fact  we  examine  Figure  7 (or  see  Table  4.2)  for  the  data  of  1 

and  P3*^3/P3‘^jL  " sees  that  when  d^/d  = 1 the  temperature  dis- 

tortion can  indeed  be  negative.  Therefore,  if  d^/d  are  choosen  between 
0.8  and  1 the  error  can  be  minimized.  However  one  must  keep  in  mind 
that  the  elimination  of  error  by  heat  capacitance  can  work  during  the 
transient  period  only, for  once  a steady  state  conduction  is  established 
the  heat  capacity  pc  will  no  longer  have  any  effect  and  overheating 
at  the  cavity  base  eventually  will  develop.  This  can  be  seen  best  from 
the  governing  equation  (1)  that  at  steady  state  the  unsteady  term  which 
contains  pc  product  is  zero  and  is  not  a parameter  affecting  the  dis- 
tortion. 
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Another  Important  fact  that  should  be  mentioned  Is  that  in  general 
the  optimum  choice  of  d^/d  ratio  for  given  ratio  does  not 

vary  very  much  with  the  variation  of  c/D  ratio.  The  Insensitivity  of 
the  optimum  d^/d  ratio  to  the  e/D  ratio  ranging  from  0.02  to  0.1  means 
that  the  distortion  of  the  temperature  is  insensitive  to  the  cavity  depth 
or  the  thickness  of  the  disk.  This  fact  was  already  pointed  out  by  Chen 
and  Danh  [11]  in  their  experiment  that  the  temperature  distortion  at  the 
base  of  the  cavity  is  more  sensitive  to  the  variation  of  the  cavity 
diameter  than  the  depth  of  the  cavity  drilled. 

As  an  example  of  a practical  application,  let  us  consider  a measure- 
ment of  the  trans:!  3nt  temperature  response  of  an  engine  block  made  of 
aluminum.  From  Figure  3b  we  know  that  aluminum  has  high  thermal  conduc- 
tivity. Therefore  copper-constantan  thermocouple  which  has  a higher 
thermal  conductivity  than  aluminum  should  be  chosen.  For  this  material 
combination  we  have  p^c^/p^^c^  »=  1. 3.  Now  if  the  thermocouple 

cavity  is  drilled  such  that  e/D  =0.1  then  from  Figure  7 interpolating 
between  “ 2 and  1 for  we  find  that  the  optimum 

d^/d  for  “ 1.69  is  approximately  0.7 

One  disadvantage  of  invoking  finite  element  analysis  is  that  the 
result  does  not  give  a clear  functional  relation  among  the  parameters 
Involved.  In  an  attempt  to  obtain  a simple  and  useful  relation  to 
relate  the  various  parameters  we  note  the  following  fact  and  result: 

(a)  the  optimum  d^/d  ratio  for  zero  temperature  distortion  is  a strong 
function  of  '^3/'^2^  and  pc  ratio  but  is  relatively  insensitive  to  the  t/D 
ratio,  (b)  from  the  theoretical  reasoning  the  d^/d  ratio  is  Independent 
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of  pc  ratio  If  the  problem  is  steady  state.  A simple  steady  one  dimen- 
sional analysis  In  which  the  thermocouple  and  th  insulation  matctial  in 
the  cavity  Is  made  to  conduct  the  same  amount  of  heat  that  would  be 
transferred  without  the  cavity  gives  the  relation 

d^/d  = - <2) 

Using  the  above  equation  as  a base  we  find  that  for  the  transient  heat 
conduction  as  calculated  by  the  finite  element  method  the  following 
equation  (10)  correlates  very  well  with  the  optimum  d^/d  ratio. 

dt/d  - (p3C3/p^c^)°-^  - <^)  (10) 

Equation  (10)  gives  an  error  or  distortion  of  no  more  than  two  per- 
centage points.  In  practice  equation  (10)  may  be  used  as  a rule  of 
thumb. 


CONCLUSION 

An  analysis  of  the  temperature  distortion  caused  by  the  cavity 
drilled  into  a disk  to  accommodate  the  thermocouple  has  been  studied. 
The  calculation  Is  carried  out  for  the  case  of  constant  heat  flux.  It 
Is  shown  that  the  temperature  at  the  base  of  the  cavity  distorted  from 
that  without  a cavity  can  be  eliminated  by  a properly  chosen  combina- 
tion of  the  ratio  of  the  thermocouple  diameter  to  the  cavity  diameter, 
d^/d  and  the  thermocouple  material  optimum  ratio  of  d^/d 

can  be  found  from  Figures  7,  8 and  9 or  Tables  2,  3 and  4,  or  approxi- 
mately from  equation  (10).  As  a rule  the  thermocouple  must  be  chosen 


15 


to  have  a higher  thermal  conductivity  than  that  of  the  heat  conducting 
solid.  The  cavity  diameter  should  be  as  small  as  practically  possible. 
For  the  case  of  time  dependent  surface  heat  flux  the  present  result  may 
be  also  used  as  a general  guide. 
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PART  I Table  la 


Connnonly  Used  Thermocouples 


Types  of  Thermocouples 

Temperature 

°F 

Ranges 

“C 

Copper /Conatantan 

-300/750 

-148/398 

Iron/Constantan 

-300/1600 

-148/871 

Chromel/Alumel 

-300/2300 

-148/1260 

Chromel/Constantan 

32/1800 

0/982 

Platinum  10%  Rhodium/Platinum 

32/2800 

0/1537 

Platinum  13%  Rhodlum/Platlnum  6% 

Rh 

100/3270 

37/1798 

Platinel  1813  Platinel  1503 

32/2372 

0/1300 

Iridium/ Iridium  60%  Rhodium 

2552/3326 

1400/1830 

Tungsten  3%  Rhenium/Tungsten  25%  Rhenium 

50/4000 

10/2204 

Tungsten/Tungsten  26%  Rhenium 

60/5072 

15/2800 

Tungsten  5%  Rhenium/Tungsten  26% 

Rhenium 

32/5000 

0/2760 
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PART  I Table  lb 


Thermal  Properties  of  Thermocouple  Materials 


cal 

"3 

'^3 

cal 

P3C3 

P3C3 

Metal,  Insulator 
(subscript  3) 

sec.  cm < 

K Steel 

K aluminum 

P — 4 

3 3 cm-* 

pc  steel  pc 

aluminum 

Aluminum 

0.554 

4.66 

1.00 

0.65 

0.60 

1.0 

Copper 

0.935 

7.87 

1.69 

0.84 

0.78 

1.30 

Chromium 

0.201 

1.70 

0.36 

0.84 

0.78 

1.30 

Nickel 

0.160 

1.35 

0.29 

1.20 

1.12 

1.86 

Platinum 

0.174 

1.46 

0.31 

0.74 

0.69 

1.14 

Steel 

0.119 

1.0 

0.21 

1.08 

1.00 

1.66 

Tungsten 

0.373 

3.14 

0.67 

0.59 

0.55 

0.92 

Iridium 

0.345 

2.91 

0.62 

0.66 

0.61 

1.01 

Rhodium 

0.360 

3.03 

0.65 

0.72 

0.67 

1.12 

Rhenium 

0.111 

0.94 

0.20 

0.73 

0.68 

1.14 

Nlc  ke 1-Chromlum 

0.040 

0.34 

0.07 

1.01 

0.94 

1.56 

Copper-Nickel 

0.059 

0.50 

0.11 

0.99 

0.92 

1.52 

Teflon 

6 X 10 

-4 

0.005 

0.001 

0.54 

0.5 

0.83 

Air 

1. 

26  X 10 

0.001 

0.0002 

0.53 

0.49 

0.82 

Note:  The  value  quoted  Is  the 

are  available. 

averaged  value  over 

200  tc  800 

°tc  whenever 

data 

PART  I Table  2 


Percentage  Error  of  Temperature  at  Cavity  Base  ( /D  = 0.1) 


Table  2.1 

1C3/K:  j 

= 0.5 

e/D  = 

0.1 

P3C3/P1C1 

0.0 

0.2 

0.4 

0.6 

0.8 

1.0 

T X 

10^ 

Error 

% 

0.9 

1 

1.9 

1.4 

1.0 

0.8 

0.6 

0.5 

2 

3.4 

2.5 

1.9 

1.4 

1.1 

1.0 

4 

5.1 

3.8 

2.9 

2.2 

1.7 

1.5 

8 

6.4 

4.7 

3.7 

2.9 

2.3 

1.9 

12 

6.9 

5.1 

4.1 

3.3 

2.6 

2.2 

20 

7.4 

5.5 

4.5 

3.7 

2.9 

2.4 

40 

7.8 

5.9 

4.8 

4.0 

3.2 

2.6 

60 

7.9 

6.0 

4.9 

4.0 

3.2 

2.6 

80 

7.9 

6.0 

4.8 

3.9 

3.2 

2.6 

100 

7.9 

6.0 

4.8 

3.9 

3.2 

2.6 

1.3 

1 

1.3 

0.9 

0.6 

0.4 

0.3 

2 

2.4 

1.7 

1.1 

0.8 

0.6 

4 

3.7 

2.6 

1.8 

1.2 

1.1 

8 

4.6 

3.4 

2.4 

1.6 

1.5 

12 

5.0 

3.8 

2.7 

1.9 

1.8 

20 

5.4 

4.2 

3.1 

2.2 

2.0 

40 

5.8 

4.5 

3.4 

2.5 

2.3 

60 

5.9 

4.5 

3.4 

2.5 

2.3 

80 

5.9 

4.5 

3.3 

2.5 

2.3 

100 

5.8 

4.5 

3.3 

2.5 

2.3 

1.8 

1 

1.2 

0.7 

0.4 

0.1 

0.1 

2 

2.3 

1.4 

0.8 

0.4 

0.3 

4 

3.6 

2.3 

1.4 

0.7 

0.6 

8 

4.5 

3.0 

1.9 

0.9 

1.0 

12 

4.9 

3.4 

2.1 

1.1 

1.2 

20 

5.3 

3.8 

2.5 

1.5 

1.6 

40 

5.6 

4.1 

2.8 

1.7 

1.9 

60 

5.7 

4.1 

2.7 

1.7 

1.9 

80 

5.7 

4.1 

2.6 

1.6 

1.9 

100 

5.7 

4.0 

2.6 

1.6 

1.9 

I 


Table  2.2 

<3/.Ci 

= 1.0 

e/D 

= 0.1 

PjCj/Pi^i 

d /d  0.0 

0.2 

0.4 

0.6 

0.8 

1.0 

T X 10^ 

Error 

% 

0.  7 

1 

1.2 

0.8 

0.5 

0.3 

0.2 

2 

2.3 

1.5 

1.0 

0.6 

0.2 

4 

3.4 

2.4 

1.6 

1.0 

0.3 

8 

4.3 

3.1 

2.1 

1.3 

0.3 

12 

4.6 

3.3 

2.3 

1.5 

0.3 

20 

4.9 

3.6 

2.6 

1.6 

0.3 

40 

5.2 

3.8 

2.8 

1.7 

0.2 

60 

5.2 

3.9 

2.8 

1.7 

0.2 

80 

5.2 

3.9 

2.8 

1.8 

0.2 

100 

5.2 

3.9 

2.8 

1.8 

0.2 

0.95 

1 

1.2 

0.7 

0.4 

0.2 

0.0 

2 

2.2 

1.4 

0.8 

0.3 

0.0 

4 

3.3 

2.2 

1.2 

0.6 

0.0 

8 

4.2 

2.8 

1.7 

0.8 

0.1 

12 

4.5 

3.1 

1.9 

0.9 

0.0 

20 

4.8 

3.4 

2.2 

1.1 

0.0 

40 

5.1 

3.6 

2.3 

1.2 

0.0 

60 

5.1 

2.6 

2.3 

1.2 

0.0 

80 

5.1 

3.6 

2.3 

1.2 

0.0 

100 

5.1 

3.6 

2.3 

1.2 

0.0 

1.3 

1 

1.1 

0.6 

0.2 

-0.0 

-0.1 

2 

2.1 

1.1 

0.5 

-0.0 

-0.2 

4 

3.2 

1.8 

0.8 

0.0 

-0.3 

8 

4.0 

2.4 

1.2 

0.2 

-0.3 

12 

4.4 

2.7 

1.4 

0.3 

-0.3 

20 

4.7 

3.0 

1.6 

0.4 

-0.3 

40 

5.0 

3.2 

1.7 

0.4 

-0.2 

60 

5.0 

3.2 

1.7 

0.4 

-0.2 

80 

5.0 

3.2 

1.6 

0.4 

-0.2 

100 

5.0 

3.1 

1.6 

0.4 

-0.2 

1 

0.3 

-0.4 

-0.8 

-1.0 

-1.4 

2 

0.7 

-0.6 

-1.4 

-1.9 

-2.7 

4 

1.1 

-0.8 

-2.0 

-2.9 

-4.6 

8 

1.3 

-0.9 

-2.4 

-3.7 

-6.8 

12 

1.3 

-0.8 

-2.3 

-3.8 

-8.3 

20 

1.2 

-0.4 

-1.8 

-3.7 

-9.9 

40 

1.4 

0.2 

-1.0 

-3.3 

-11.1 

60 

1.4 

0.3 

-0.9 

-3.2 

-11.3 

80 

1.4 

0.4 

-0.8 

-3.2 

-11.3 

100 

1.4 

0.4 

-0.8 

-3.2 

-11.3 

\ \ ' 

! . • 
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PART  I Table  3 


Percentage  Error  of  Temperature  at  Cavity  ( /D  = 0.04) 


Table  3.1 

K3/K1 

= 0.5 

e/D  = 

0.04 

P3C3/P1C1 

d^/d 

0.0 

0.2 

0.4 

0.6 

0.8 

1.0 

T X 10 

2 

Error 

% 

0.9 

.2 

1.5 

0.9 

0,7 

0.5 

0.4 

0.4 

.4 

3.2 

2.0 

1.5 

1.1 

0,8 

0.7 

.8 

5.6 

3.9 

2.8 

2,0 

1.5 

1.3 

1.6 

8.2 

6.0 

4.5 

3.3 

2.4 

1.9 

2.4 

9.3 

7.0 

5.3 

3,9 

2.9 

2.3 

4.0 

10.3 

7.8 

6.1 

4.5 

3.4 

2.7 

8.0 

11.2 

8.5 

6.7 

5.2 

3.9 

3.1 

12.0 

11.5 

8.8 

7.1 

5.5 

4.2 

3.3 

16.0 

11.8 

9.0 

7.3 

5.8 

4.5 

3.4 

20.0 

11.9 

9.1 

7.5 

6.0 

4,6 

3.5 

1.3 

.2 

0.8 

0.5 

0.3 

0.2 

0.2 

.4 

1.9 

1,2 

0,8 

0.5 

0.4 

.8 

3.7 

2.5 

1.6 

1.0 

0.8 

1.6 

5,8 

4.0 

2.6 

1.6 

1.3 

2.4 

6.8 

4.9 

3.2 

2.1 

1.6 

4.0 

7.6 

5.6 

3.8 

2,5 

2.0 

8.0 

8.4 

6.3 

4.4 

2.9 

2.5 

12.0 

8.6 

6.6 

4.7 

3.2 

2.7 

16.0 

8.8 

6.9 

5.0 

3.4 

2,9 

20.0 

9.0 

7.0 

5.2 

3.6 

3.0 

1.8 

.2 

0.7 

0.4 

0.2 

0,1 

0.0 

.4 

1.7 

0.9 

0.5 

0,2 

0,1 

.8 

3.4 

2.1 

1.1 

0,5 

0.3 

1.6 

5.5 

3.5 

2.0 

0.9 

0.6 

2.4 

6.5 

4.3 

2.5 

1.2 

0.9 

4.0 

7.4 

5.1 

3.1 

1.6 

1.2 

8.0 

8.2 

5.8 

3.6 

1.9 

1.7 

12.0 

8.5 

6.1 

3.9 

2,1 

2.0 

16.0 

8.7 

6.3 

4.1 

2.3 

2.2 

20.0 

8.8 

6.5 

4.3 

2.5 

2,4 

4..~ 


Table  1.2 


c/D  = 0,04 


d /d  0.0 

t 

0.2 

0.4 

0.6 

0.8 

1.0 

r 

T X 10^ 

Error 

% 

1 0.7 

0.2 

0.8 

0.5 

0.3 

0.2 

0.1 

1 

0.4 

1.8 

1.1 

0.7 

0.4 

0.3 

0.8 

3.4 

2.2 

1.4 

0.8 

0.4 

1.6 

5.4 

3.6 

2.2 

1.3 

0.4 

• 

2.4 

6.3 

4.  3 

2.7 

1,6 

0.4 

4.0 

7.1 

5.0 

3.2 

1.9 

0.4 

8.0 

7.6 

5.5 

3.7 

2.2 

0.4 

12.0 

7.9 

5.8 

4.0 

2.4 

0.4 

16.0 

8.0 

5.9 

4,1 

2.4 

0.4 

[ 

20,0 

8.1 

6.1 

4.2 

2,5 

0,3 

1 

[ 0.95 

0.2 

0.7 

0.4 

0.2 

0.1 

0.0 

1 ' 

0.4 

1.6 

0.9 

0.5 

0.2 

0,0 

1 ‘ 

0.8 

3.3 

1.9 

1.0 

0.4 

0.1 

f ' ' 

1.6 

5.2 

3,3 

1,8 

0,7 

0.1 

I 

2.4 

6.1 

3.9 

2,2 

0.9 

0.1 

4.0 

6.9 

4.6 

2.6 

1.1 

0.1 

} 

8.0 

7.5 

5,1 

3.1 

1.4 

0,1 

12.0 

7.  7 

5,4 

3.3 

1,6 

0.1 

16.0 

7.9 

5.6 

3.5 

1,6 

0,1 

20.0 

8.0 

5,7 

3.6 

1.7 

0.1 

t ’ 1-3 

0,2 

0.6 

0.3 

0.0 

-0.1 

-0.1 

? * 

0.4 

1.5 

0.7 

0.2 

-0.1 

-0.2 

' ■ 1 

0.8 

3.1 

1,6 

0.6 

-0.1 

-0,3 

! 

1.6 

5.0 

2.8 

1.2 

0.1 

-0.4 

i ■ 

2.4 

5.9 

3.5 

1.5 

0.2 

-0.4 

' 

1 

4.0 

6.7 

4.1 

1.0 

0.3 

-0.4 

- i 

8.0 

7.3 

4.6 

2.2 

0.4 

-0.4 

f 

12.0 

7.6 

4.9 

2.5 

0.5 

-0.4 

t 

16.0 

7.7 

5.1 

2.7 

0.6 

-0.4 

20.0 

7.9 

5.2 

2,8 

0.7 

-0.3 

r 


Table  3.3 

K3/K1 

= 5 

e/D 

= 0.04 

P3C3/P1C1 

d^/d  0.0 

0.2 

0.4 

0.6 

0.8 

.10 

T X 10^ 

Error 

% 

0.6 

0.2 

0.3 

-0.1 

-0.3 

-0.4 

-0.6 

0.4 

0.9 

0.0 

-0.5 

-0.8 

-1.1 

0.8 

2.0 

0.4 

-0.6 

-1.2 

-2.1 

1.6 

3.4 

0.9 

-0.7 

-1.8 

-3.2 

2.4 

4.1 

1.2 

-0.8 

-2.1 

-4.0 

4.0 

4.6 

1.4 

-0.9 

-2.5 

-5.0 

8.0 

4.7 

1.4 

-1.1 

-3.0 

-6.4 

12.0 

4.7 

1.4 

-1.0 

-3.2 

-7.3 

16.0 

4.7 

1.6 

-0.8 

-3.2 

-7.9 

20.0 

4.7 

1.7 

-0.6 

-3.2 

-8.3 

Table  3.4 

K3/K1  = 10 

e/D 

= 0.04 

0.75 

0.2 

0.0 

-0.4 

-0.6 

-0.7 

-0.9 

0.4 

0.3 

-0.6 

-1.1 

-1.4 

-1.7 

0.8 

1.0 

-0.7 

-1.7 

-2.3 

-3.0 

1.6 

2.1 

-0.8 

-2.5 

-3.5 

-4.6 

2.4 

2.6 

-0.8 

-2.9 

-4.2 

-5.7 

4.0 

2.9 

-0.9 

-3.4 

-5.0 

-7.2 

8.0 

2.9 

-1.0 

-3.7 

-5.8 

-9.5 

12.0 

2.8 

-0.8 

-3.4 

-5.8 

-11.0 

16.0 

2.8 

-0.5 

-3.1 

-5.7 

-11.9 

20.0 

2.8 

-0.3 

-2.7 

-5.6 

-12.5 

-'•'V . 


PART  I Table  A 


Percentage  Error  of  Temperature  at  Cavity  Base  (c/D  = 0.02) 


Table  4.1 

K^l  < 

1 = 

e/D  - 

0.02 

P3C3/P1C1 

AJA  0.0 

0.2 

0.4 

0.6 

0.8 

1.0 

T X 10 

Error 

% 

1.3 

0,05 

0.3 

0.1 

0.1 

0.1 

0.0 

0.10 

0.8 

0.4 

0,2 

0.2 

0.1 

0.20 

2.0 

1.2 

0,7 

0.5 

0,4 

0.40 

4,5 

2.7 

1,6 

1.0 

0.8 

0.60 

6,5 

4,0 

2.  3 

0.4 

1,1 

1.00 

8,9 

5.6 

3.3 

2.0 

1.5 

2.00 

11.4 

7.6 

4.7 

2.8 

2,2 

3.00 

12.4 

8.5 

5.4 

3.3 

2.6 

4.00 

12.9 

9.0 

5.8 

3.6 

2.8 

5,00 

13.1 

9.3 

6.0 

3.7 

3.0 
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Table  A.  2 

k^/k 

1 = ' 

c/D 

= 0.02 

PsC^/PlC^ 

d^/d  0.0 

0.2 

0.4 

0.6 

0,8 

1.0 

T X 10^ 

Error  °l 

0.7 

0.05 

0.3 

0.2 

0.1 

0.1 

0.1 

0.10 

0.8 

0.5 

0,3 

0,2 

0.2 

0.20 

1.9 

1.1 

0,7 

0.4 

0,3 

0.40 

4.3 

2.5 

1.5 

0,8 

0.5 

0.60 

6.1 

3.6 

2.1 

1.1 

0.5 

1.00 

8.3 

5.1 

2.9 

1.5 

0.6 

2.00 

10.6 

6.8 

3.9 

2,1 

0.6 

3.00 

11.4 

7.5 

4.4 

2.4 

0.5 

* 

4.00 

11.8 

7.8 

4.7 

2,5 

0,5 

5.00 

12.1 

8.1 

4.9 

2.8 

0.5 

0.95 

0.05 

0.2 

0.1 

0.0 

0.0 

0.0 

0.10 

0.7 

0.3 

0.1 

0.1 

0.0 

0.20 

1.8 

0.9 

0.4 

0.1 

0.0 

0.40 

4.0 

2.1 

1.0 

0.3 

0.1 

0.60 

5.8 

3.1 

1.5 

0.5 

0.1 

, 

1.00 

8.0 

4.5 

2.2 

0.8 

0.1 

] 

2.00 

10.3 

6.1 

3,1 

1,2 

0.1 

. j 

3.00 

11.2 

6.8 

3.5 

1.4 

0.1 

4.00 

11.6 

7.2 

3.8 

1.5 

0.1 

i 

5.00 

11.8 

7.4 

4.0 

1.7 

0.1 

' 1.3 

0.05 

0.2 

0.0 

-0.0 

-0.1 

-0.1 

0.10 

0.5 

0.1 

-0.1 

-0.1 

-0.1 

* 

0.20 

1.6 

0.6 

0.1 

-0.2 

-0,3 

i 

0.40 

3.7 

1.6 

0.5 

-0.2 

-0.4 

V ' 

0.60 

5.4 

2,5 

0.8 

-0.1 

-0.5 

fe 

1.00 

7.6 

3.8 

1.4 

-0.0 

-0.5 

1 

2.00 

10.0 

5.4 

2.1 

0.2 

-0.5 

1 

1 

3.00 

10.9 

6.0 

2.5 

0,3 

-0.5 

1 

1 

4.00 

11.3 

6.4 

2.7 

0.4 

-0.5 

1 

1 

5.00 

11.5 

6,7 

2.9 

0.4 

-0.5 

Table  4.3 

= 10 

e/D 

= 0.02 

P3C3/P1C1 

AjA 

0 

0.2 

0.4 

0.6 

0.8 

1.0 

^ X 10 

2 

Error  % 

0.75 

0.05 

0.7 

0.0 

-0.26 

-0.33 

-0.37 

-0.41 

0.10 

1.7 

0.0 

-0.48 

-0.66 

-0.76 

-0.86 

0.20 

4.9 

0.2 

-0.72 

-0.93 

-1.40 

-1.61 

0.40 

9.0 

1.2 

-0.85 

-1.78 

-2.28 

-2.71 

0.60 

11.6 

2.0 

-0.86 

-2.20 

-2.89 

-2.50 

1.00 

13.2 

3.2 

-0.84 

-2.78 

-3.76 

-4.66 

2.00 

15.8 

4.5 

-0.88 

-3.66 

-5.06 

-6.44 

3.00 

16.7 

4.9 

-0.99 

-4.20 

-5.84 

-7.56 

4.00 

17.1 

5.1 

-1.10 

-4.57 

-6.38 

-8.42 

5.00 

17.4 

5.1 

-1.20 

-4.80 

-6.75 

-9.12 

Table  4.4 

K3/K1 

= 5.0 

c/D  = 

0.02 

0.6 

0.05 

0.7 

0.05 

-0.00 

-0.17 

-0.20 

-0.25 

! 

0.10 

1.7 

0.24 

-0.15 

-0.32 

-0.43 

-0.55 

1 

0.20 

4.9 

0.87 

-0.06 

-0.53 

-0.80 

-1.09 

0.40 

9.0 

1.33 

0.34 

-0.70 

-1.29 

-1.89 

1 

0.60 

11.6 

3.53 

0.71 

-0.78 

-1.61 

-2.48 

1.00 

13.2 

5.12 

1.25 

-0.88 

-2.04 

-3.33 

1 

2.00 

15.8 

6.87 

1.92 

-1.00 

-2.62 

-4.59 

3.00 

16.7 

7.48 

2.15 

-1.12 

-2.97 

-5.33 

4.00 

17.1 

7.72 

2.21 

-1.25 

-3.26 

-5.88 

5.00 

17.4 

7.80 

2.20 

-1.36 

-3.49 

-6.32 

V ' 
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PART  II 


IMPROVED  ACCURACY  IN  THE  PREDICTION  OF  SURFACE  IffiAT  FLUX 
AND  TEMPERATURF  BY  AN  INTRINSIC  THERMOCOUPLE 
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INTRODUCTION 


In  the  study  of  transient  heat  transfer,  many  experimental  difficulties 
may  arise  If  heat  flux  sensors  or  thermocouples  are  Installed  direct  at  the 
surface  of  a body.  For  example,  a probe  may  be  damaged  by  a piston  or  a 
projectile  sliding  over  a cylinder  or  barrel.  A probe  on  a melting  and 
ablative  surface  of  heat  shield  can  be  easily  destroyed  because  of  high 
temperature.  Furthermore  a surface  probe  exposed  to  both  radiative  and 
convective  environment  may  measure  an  erroneous  surface  heat  flux  and 
temperature  If  the  probe  has  a different  radiative  property  from  that  of 
the  measured  surface.  In  these  circumstances,  calculation  of  the  transient 
surface  heat  flux  and  the  surface  temperature  can  be  achieved  by  Inverting 
a temperature  history  measured  at  some  location  Inside  the  body. 

In  general,  the  prediction  of  a surface  heat  flux  and  temperature  by 
the  measured  data  at  some  location  Interior  to  a body  Is  known  as  the 
"Inverse  problem".  Many  configurations,  such  as  spheres,  cylinders,  and 
slabs,  had  been  studied  by  many  workers  and  many  methods  such  as  numerical, 
graphical,  series,  convolution  Integral,  and  Laplace  transforms  were  used. 
Stolz  [1],  Beck  [2]  and  Williams  and  Curry  [3],  considered  the  numerical 
Inversion  of  the  Integral  solution  for  seml-lnf Inlte  and  other  bodies. 

In  this  method,  care  Is  required  In  selecting  a time  Interval  In  order 
to  achieve  a stable  solution.  Carslaw  and  Jaeger  [A],  Burggraf  [5], 
Koveryanov  [6],  and  Shumakov  [7],  respectively  considered  different  series 
approaches  In  which  generally  the  local  heat  flux  at  an  Interior  location 
and  their  higher  derlvatles  are  required.  However,  It  Is  difficult  to 
measure  experimentally  or  to  process  the  measured  data  for  the  derivative 
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of  the  temperature.  Sparrow,  Hajl-Sheikh,  and  Lundgren  [8],  Imber  and 
Kahn  [9],  Imber  [10],  Sabherwal  [11],  Masket  and  Vastano  [12],  Deverall  and 
Channapragada  [13]  and  Chen  and  Thomsen  [14]  applied  the  transform  method. 

In  these  works,  the  solution  Is  represented  In  either  an  integral  form 
after  some  manipulation  of  the  contour  integral  from  the  Inverse  trans- 
form, or  in  a series  form  after  an  expansion  of  the  solution  for  small  and 
large  times.  Using  Laplace  transformation  Chen  and  Thomsen  [14]  introduced 
a polynomial  in  terms  of  an  error  function  to  represent  the  response  of 
thermocouple  measurement  and  the  inversion  is  accomplished  for  any  transient 
surface  heat  flux  at  the  inner  surface  of  a cylindrical  tube.  In  their 
study,  the  cylindrical  thickness  was  assumed  to  be  relatively  thick  such 
that  the  temperature  at  a large  distance  from  the  heating  surface  remains 
constant.  Therefore,  only  one  interior  temperature  response  near  the  sur- 
face was  needed  in  the  experimental  measurement.  Their  inversion  solution 
however  was  valid  only  for  a short  duration  due  to  the  asymptotic  expansion 
of  the  modified  Bessel  function  in  the  inverse  Laplace  transform.  Chen  and 
Chlou  [15]  studied  the  inversion  problem  for  the  case  of  a semi-infinite 
slab  or  a thick  slab  using  a Laplace  transformation.  The  exact  solution 
was  obtained  from  the  inverse  Laplace  transform  for  any  time  interval.  It 
was  then  shown  that  their  analysis  may  be  approximately  applied  to  the  case 
of  the  hollow  cylinder  if  the  interior  temperature  response  is  measured  at 
a location  close  to  the  inner  wall. 

This  report  presents  (a)  the  improved  numerical  solution  of  the  in- 
version solution  reported  by  Chen  and  Chiou  [15]  and  (b)  a further  demon- 
stration of  the  capability  of  the  solution. 
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The  theoretical  analysis  of  Chen  and  Chiou  [15]  Is  recapitulated  in 
Appendix  A In  which  the  surface  heat  flux  and  temperature  is  predicted  by 
inverting  a temperature  history  measured  at  some  location  inside  the  solid 
body.  The  inversion  solution  is  obtained  by  invoking  Laplace  transformation. 
Both  the  surface  heat  flux  and  temperature  are  given  by  Eqs.  (19)  and  (20) 
in  Appendix  A. 

It  was  thought  that  the  accuracy  of  the  computer  program  generated 
for  the  solution  in  the  previous  report  by  Chen  and  Chiou  [15]  can  be  im- 
proved further  for  the  following  reasons.  First,  the  coefficients  b 

n 

(see  Appendix  A Eq.(ll))  in  the  previous  formulation  has  a dimension  of 

temperature.  Therefore  the  determination  of  the  coefficients  depends  on 

the  temperature  range  of  each  particular  experiment.  It  was  found  tliat  the 

absolute  value  of  the  coefficients  b in  some  cases  can  become  as  large  as 

n 
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an  order  of  10  . Therefore  during  the  subsequently  numerical  manipulation 

in  the  computer  program  error  due  to  round  off  and  the  standard  fixed  up 
when  an  underflow  occurred  may  become  appreciable.  To  remedy  this  diffi- 
culty the  dimensionless  formulation  is  introduced  in  the  analysis  (Appendix 

A)  in  which  the  coefficient  b is  also  made  dimensionless.  As  a result 

n 

the  magnitude  of  the  coefficient  b^  can  be  greatly  reduced.  Secondly, 
the  double  precision  format  was  not  used  throughout  the  previous  computer 
program.  It  is  felt  that  further  accurate  results  may  be  obtained  if  the 
double  precision  format  is  adopted  in  the  program. 

In  the  following  section  the  new  computer  solution  is  shown  to  be 
indeed  more  accurate.  Later  the  solution  is  shown  to  be  capable  of  pre- 
dicting a case  involving  a pei iodic  surface  heat  flux  or  periodic  tempera- 
ture variation. 
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RESULTS  OF  THE  IMPROVED  COMPUTER  PROGRAM 


The  previous  computer  program  of  Chen  and  Chlou  [15]  was  recast  in 
dimensionless  form  and  written  In  the  double  precision  format.  The  new 
computer  program  is  listed  In  Appendix  B.  The  results  predicted  by  the 
new  and  previous  computer  program  are  given  in  Appendix  C and  shown  in 
Figure  1 for  the  case  of  the  constant  surface  heat  flux.  This  is  the  case 
in  which  a steel  slab  initially  at  a uniform  temperature  is  suddenly  sub- 
jected to  a constant  heat  flux  Q at  one  of  the  surfaces  and  kept  at  the 
initial  temperature  on  the  other  surface.  Figure  1 shows  the  solution 
predicted  by  inverting  the  temperature  response  at  an  interior  of  the  slab 
from  the  new  and  previous  computer  program.  This  solution  predicted  by 
the  new  and  previous  programs  used  the  ten  term  representation  for  the 
thermocouple  response.  The  comparison  clearly  shows  the  improvement  of 
the  new  solution  over  the  previous  one.  Except  for  the  short  time  duration 
the  solution  with  the  new  program  reduces  the  error  to  only  one  half  of 

the  error  of  the  previous  program  i.e.,  an  error  of  less  than  one  percent. 

* 

In  the  short  time  period  the  solution  exhibits  a Gibbs  phenomenon  because 

of  the  discontinuity  of  the  surface  temperature  gradient  occurred  at  initial 

condition.  The  solution  shows  a 17%  of  initial  overshot  of  heat  flux  and 

then  a 7.8%  of  undershot  before  the  solution  approaches  the  constant  heat 

flux.  It  should  be  remarked  that  Gibbs  phenomenon  is  artificially 
★ 

Gibbs  phenomena  [3]:  for  a sequence  of  transformation  Tn(t),  n = 1,  2,... 

of  a function  q(t)  (here  q(t)  = constant)  if  the  interval  Tn(t), 

n*^oo 

Tn(t)]  contains  points  outside  the  interval  inf  q(t), 

t*^c  ri'^oo 

. . o o 

11m 

t-*-t^  sup  q(t)]  then  the  sequence  is  said  to  exhibit  a Gibbs  phenomena. 
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introduced  due  to  the  idealization  of  the  initial  condition.  In  most 
|)ractli:al  Hlluations  the  surface  l»eut  flux  will  bo  continuous.  Thorol  on* , 
Ibbs  phenomenon  will  not  appear. 

Figure  2 (see  Table  1 also)  shows  the  comparison  between  tlie  solutions 
for  the  constant  heat  flux  case  with  10  and  20  term  representation  for  the 
thermocouple  response.  One  sees  that  the  solution  with  20  term  representa- 
tion after  the  initial  Gibbs  phenomenon  quickly  approaches  the  expected 
constant  heat  flux  solution  with  a negligible  error  of  less  than  0.14  per- 
cent. This  shows  the  accuracy  of  the  new  computer  program.  From  Figure  2 
one  also  observes  that  both  overshoot  and  undershoot  of  Gibbs  phenomenon 
are  smaller  for  the  20  term  representation.  Additionally  the  points  of  the 
overshoot  and  the  undershoot  have  moved  to  near  the  zero  time  which  agrees 
with  the  characteristic  of  Gibbs  phenomenon.  According  to  Gibbs  phenomenon 
the  point  of  overshoot  should  approach  the  initial  zero  if  the  number 
of  terms  of  the  series  which  represent  the  thermocouple  response  is  increased 
to  infinite. 


VERIFICATION  OF  OSCILLATORY  SOLUTION 

As  a measure  of  applicability  of  the  present  inversion  solution,  a 
test  problem  was  solved  for  the  case  of  a slab  subjected  to  a periodic 
surface  temperature  variation  on  one  surface  and  held  to  the  initial 
temperature  on  the  other.  The  analytic  solution  for  the  problem  is  given 
in  Appendix  IV  where  a more  suitable  form  of  the  solution  than  the  one 
given  by  Carslaw  and  Jaeger  [4]  is  derived  and  tabulated  for  the  thermo- 
couple response  at  one  tenth  of  the  slab  thickness  from  the  surface.  The 


PART  II  TAP.I.F  1 

Comparlsiou  of  Inversion  Prediction  and  Exact  Solution 


0(1, t 
0.1  0.0624 

0.2  0.1628 
0.4  0.3395 

0.6  0.4882 

0.8  0.6183 

1.0  0.7353 


0.1 

0.0624 

0.2 

0.1628 

0.4 

0.3395 

0.6 

0.4882 

0.8 

0.6183 

1.0 

0.7353 

6(o,t)~r  e(o,t)  TMror-^ 

At 

EXACT  PREDICTED  % 

0.3758  0.3773  +0.399 

0.5315  0.5314  +0.019 

0.7516  0.7516  0.00 

0.9205  0.9205  0.00 

1.0629  1.0628  -0.009 

1.1884  1.1878  -0.05 

0.3758  0.4048  +7.71 

0.5315  0.5284  -0.58 

0.7516  0.7516  0.00 

0.9205  0.9205  0.00 

1.0629  1.0626  -0.028 

1.1884  1.1887  +0.025 


30  (0,t) 
3x 

EXACT 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 


^0  (0,t)  ERROR 
3x 

PREDICTED  X 


0.9873  -1.27 
1.0014  +0.14 
1.0001  +0.01 
1.0000  0.00 
0.9996  -0.04 
1.9986  -0.14 


1.0000 

1.0796 

1.0000 

0.9391 

1 . 0000 

1.0179 

1.0000 

0.9923 

1.0000 

1,0043 

1.0000 

0.9963 

ERROR  t * ((PREDICTED) -(EXACT)) /(EXACT) 


siirfart'  Is  sub  [I'ctfd  to  a perloilic  lomporatiiro  variation  wltii  a porliui  of 


H ml  I I 1 soi’t'oiHis.  l•'Ifll'on  data  points  of  the  temiuTat  arc  rt-sponso  al  l ho 
thermocouple  location  are  then  Input  to  the  inversion  pro(;ram  lOr  prediction 
of  tile  surface  temperature  and  heat  flux.  The  result  is  shown  in  I'igure  1 
and  Appendix  C where,  the  data  symbol  "2"  denotes  the  thermocouple  re- 
sponse and  "1"  the  surface  temperature.  The  accuracy  of  the  inversion  pro- 
gram is  shown  in  Table  2.  Except  for  the  extremely  short  time  period  of 
0.4  milliseconds  the  prediction  by  the  inversion  program  with  15  term  repn-- 
sentation  is  within  2 percent  of  error. 

The  accuracy  can  be  improved  more  if  more  data  points  are  used. 

Figure  4 shows  the  predicted  surface  heat  flux  which  we  were  unable  to 
compute  from  the  series  solution  (Eq.  (7)  of  Appendix  D) . This  demonstrates 
the  versatility  of  the  inversion  solution. 


APPLICATION  OF  THE  INVERSION  PROGRAM 


Three  sets  of  data  (see  Appendix  C)  provided  by  Rock  Island 
Arsenal  for  the  temnerature  response  of  a thermocouple  embedded  in  a M60 
gun  barrel  were  utilized  to  evaluate  the  inversion  solution.  The  inversion 
prediction  for  the  surface  heat  flux  from  all  three  sets  of  data  were  ex- 
tremely high  when  compared  with  other  known  data  calculated  by  Chen  and 
Chiou  [15].  Since  the  program  correctly  predicted  the  surface  heat  flux 
for  other  sets  of  experimental  data  it  was  judged  that  the  three  sets  of 
data  may  contain  inaccurate  initial  time.  For  most  experimentations  tin- 
recording  instrument  is  likely  to  experience  some  delay  in  responding  to 
the  extremely  fast  transient  heat  flux  typical  in  gun  bores.  Therefore 


an  advanced  shift  of  time  of  ! milliseconds  in  the  data  was  tested.  The 
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Table  2 


Comparison  of  Inversion  Prediction  and  Exact  Solution 


Time 

Surface 

Surface 

Error* 

(sec. ) 

Temperature 

Temperature 

Z 

(theoretical) 

(predicted) 

0.0002 

189.5041 

193.4060 

+3.563 

0.0004 

296.3119 

301.4311 

+2.367 

0.0006 

397.7934 

403.8998 

+1.921 

0.0008 

491.4497 

498.3988 

+1.689 

0.0010 

574.9747 

582.5683 

+1.534 

0.0012 

646.3119 

654.3358 

+1.417 

0.0014 

703.7046 

711.9396 

+1.321 

0.0016 

745.7396 

753.9650 

+1.236 

0.0018 

771.3818 

779.3802 

+1.157 

0.0020 

780.0000 

787.5615 

+1.08 

0.0022 

771.3818 

778.3093 

+1.00 

0.0024 

745.7396 

751.8527 

+0.918 

0.0026 

703.7046 

408.8443 

+0.874 

0.0028 

646.3119 

650.3441 

+0.712 

0.0030 

574.9747 

577.7933 

+0.569 

0.0032 

491.6697 

492.9790 

+0.372 

0.0034 

397.7934 

397.9893 

+0.0616 

0.0036 

296.3119 

295.2761 

-0.679 

0.0038 

189.5041 

188.8522 

-0.595 

*Error  = 

((Predicted  - 

(Exact))/ (Exact)) x 

100 
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result  Is  shown  In  Appendix  C-3.  Figures  5,  6 and  7 respectively  show 

2 

the  predicted  maximum  heat  flux  of  1650,  1235  and  2695  Btu/ft  sec.  approx- 
imately at  2 milliseconds  after  the  firing,  while  Figures  8,  9 and  10  show 
that  the  surface  temperatures  are  273,  234  and  396  degree  F . approximately 
at  6 milliseconds  after  the  firing. 
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CONCLUSIONS  AND  SUGGESTIONS 


The  new  inversion  computer  program  was  thoroughly  tested  for  its 
applicability  and  accuracy.  The  program  can  invert  an  intrinsic  tem- 
perature response  to  predict  the  case  of  a constant  surface  heat  flux 
or  the  case  of  a periodic  surface  temperature  within  2%  of  deviation 
except  at  the  extremely  short  time  period.  This  was  achieved  with  the 
maximum  number  of  20  input  data  points. 

It  is  expected  that  the  accuracy  of  the  prediction  will  Increase 
if  the  number  of  data  input  is  increased. 

Based  on  the  experience  gained  in  working  with  the  program  the 
following  suggestions  are  thought  relevant. 

(1)  In  conducting  an  experiment  it  is  vital  that  both  the  tempera- 
ture and  the  time  in  the  intrinsic  measurement  of  surface 
heat  flux  and  temperature  must  be  much  more  accurate  than 
the  direct  measurement.  This  is  because  the  inversion 
problem  always  Involves  predicting  a large  heat  flux  or 
temperature  variation  from  the  data  with  small  variation. 
Therefore , for  a slight  error  in  the  time  or  temperature 
measurement  a large  error  will  result  in  the  prediction 
of  the  surface  heat  flux  or  temperature.  For  example,  in 
the  prediction  of  gun  barrel  heat  flux  considered  in 
an  error  of  two  milliseconds  in  time  will  lead  to  100  per- 
cent error  in  prediction  of  the  surface  heat  flux. 
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(2)  In  selecting  the  data  points  for  Input  to  the  computer  program 
care  must  be  exercised  not  to  create  a locally  abrupt  jump  In  the  data 
An  abrupt  change  In  data  points  will  often  Introduce  an  abnormal  fitting 
of  a curve  In  Its  neighborhood  and  hence  resulting  In  an  Incorrect  pre- 
diction of  the  surface  heat  flux  and  temperature.  If  indeed  the  abrupt 
jump  of  the  data  must  be  used,  then  more  data  points  in  its  neighborhood 
must  also  be  chosen. 
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APPENDIX  A ANALYSIS  OF  THE  INVERSION  PROBLEM 


Consider  a slab,  having  a sufficient  wall  thickness,  L,  such  that  the 
outer  surface  temperature  has  a negligible  response  when  the  inner  surface 
is  exposed  to  a transient  heat  flux-  A probe,  for  example  a thermocouple, 
is  located  at  X = and  it  is  normally  desirable,  as  reported  by  Chen  and 
Li  [16],  to  be  close  to  the  heating  surface  since  a better  transient 
response  and  more  accurate  experimental  measurements  can  be  obtained 
to  reduce  error  amplification  in  the  mathematical  inversion  program. 

Under  these  circumstances  we  thus  assume  in  the  analysis  L/X^^  >>  1.  The 
governing  equation  for  the  transient  heat  conduction  may  be  written  in  a 
dimensionless  form  as 
2 

9 9 9 9 

^ 2 
9t  3x 

with  the  Initial  and  boundary  conditions 


e(x,  o)  “ 0 

(2) 

e(“>,  t)  = 0 

(3) 

9(1,  t)  = f(t) 

(4) 

Where  x = X/X^  is  a dimensionless  distance  from  the  heating  surface  and 
t = ot/X^  is  a dimensionless  time  or  Fourier  number  with  a being  the 
thermal  dlffuslvity  and  t the  real  time.  9 = (T  - T^)/T^  is  the  dimen- 
sionless temperature  above  the  Initially  uniform  temperature  T^  . 
f(t)  = (F(t)  - T )/T  is  the  dimensionless  measured  temperature  response 
at  X = 1 with  F(t)  being  the  measured  absolute  temperature.  The 
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inversion  problem  Is  then  given  the  interior  temperature  f(t)  to  predict 

the  surface  temperature  0(0,t)  and  the  surface  heat  flux  per  unit  area 

q(0,t)  or  — = -q(0,t)X, /T  k . Here  < is  the  thermal  conductivity 

ox  1 o 

o 

of  the  solid. 


The  above  problem  may  be  solved  by  Laplace  Transformation.  Let  the 
transformation  be: 


Kx,s)  - r 0(x,t) 


-ts 

e dt 


where  0 is  continuous  otherwise  satisfied  the  Dirlchlet's  condition.  The 
temperature  function  9 is  recovered  by  inversion  of  the  Laplace  Trans- 
formation as: 


0(x,t) 


St  , 
9e  ds 


where  c is  a suitable  positive  value.  Equation  (1)  and  the  initial 
condition  (2)  under  transformation  (5)  become: 


which  has  a solution,  with  integration  constants  A and  B, 


T/  \ A i „ -t^x 

0(x,s)  ■ A e + Be 


The  transformation  of  boundary  conditions  (3)  and  (4)  into  the  Laplace 


plane  give  0(",s)  = 0,  0(1, s)  = f(s).  Substitution  of  the  boundary 
conditions  into  equation  (8),  we  get  9 and  its  derivative  as: 


6(x,s)  = f(s)  e 


<^(1  - x) 


^ f (s)  e*^^^  - 


According  to  Chen  and  Thomsen  [14],  we  choose  the  temperature  response  as 
measured  by  a probe  to  be  represented  by  a polynomial 


f(t)  “2  K r(n+l)  i^"  erfc  (— ) 


or  in  Laplace  plane 


T(S)  . e-^  £ r<“  » -tJ- 


The  are  coefficients  of  the  expansion  to  be  determined  such  that  the 

N term  polynomial  describes  the  temperature  response  f(t)  measured  at 
X = 1.  With  equation  (12),  equations  (9)  and  (10)  can  be  simplified  to 


9(x,s)  = e 


. 2 r(n  + 


1/2+n 


The  inversion  of  equations  (13)  and  (14)  at  x = 0 give: 


.f- 


N 

e(o,t)  = ^ h t 

n=l 


n 


(15) 


- 39 (x.t) 


3x 


.Y  b r(n-H) 
X.0  h " 


(16) 


9 0 (0  t ) 

where  0(0, t)  gives  the  surface  temperature  and gives  the  surface 

q^i 

heat  flux  —=r~  as  a function  of  time. 

kT 

o 

The  integral  of  the  error  function  in  (11)  is  defined  as: 

erfc  (y)  dy 


i^”  erfc(- 


(-^)  » f 1^""^ 

2rt 

2/t 


(17) 


with  n = 0,  erfc  (y) 


f 

/IT 


-X 


e dx  . r(n)  in  equations  (11)  and  (16) 


is  the  gamma  function  or  Euler's  integral  function  of  the  second  kind. 


r(n)  = 


= S" 


-u)  n-1  , 
e 01  du 


(18) 


It  should  be  remarked  that  the  choice  of  the  particular  form  (11)  is  to 

ensure  the  convergence  of  the  solution  on  the  Laplace  plane  and  an  analytic 

inversion  back  to  the  physical  plane.  With  b^  coefficients  determined 

from  equation  (11)  and  the  experimental  measurement  of  the  temperature 

response  f(t)  at  x = 1,  the  surface  temperature,  T (t),  is  obtained 

w 

from  equation  (15)  as 


N 

r (t)  = T I 1 + ^ b t") 


(19) 
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from  equation  (16)  as: 


and  the  heat  flux,  q(t)  = 


-T  K DO 
o 

3x 


x=0 


-T  K 


q(t)  = 


^ E 


1 n=l 


b t 
n 


n-1/2  r(n+l) 
r(n+l/2) 


(20) 


The  above  solution  is  the  exact  solution  for  predicting  the  transient 
surface  heat  flux  and  temperature  valid  for  both  short  and  long  time  dura- 
tions as  long  as  the  slab  is  thick  enough  such  that  the  outer  surface 
■ „ maintains  its  initial  temperature.  The  feature  of  the  present  solution  is 

the  polynomial  (11)  which  on  the  Laplace  plane  gives  a term  in  equation 
I “ (12)  exp  (-/s)  to  cancel  the  term  exp  (/s)  in  equations  (9)  and  (10). 


This  polynomial  (11)  as  suggested  by  Chen  and  Thomsen  [15]  makes  the  pres- 
ent solution  simpler  than  many  inversion  solutions  derived  in  the  past  and 
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THE  SLAB  MlTH  PERICOIC  SUKEACE  TCHPEBATURE 

temperature  and  heat  ELUA 

PRECICTEO  BY  HSAI-YIN  LEE 
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THE  SLA9  tflTH  PERIOCIC  SUREACE  TEMPERATURE 


APPENDIX  C - 3 


(OMLf SSJ 

C.03os;>7 
3.073SI'* 
C.ll0£72 
C.  U7£2S 
C.  le^V  12b 
C.2?l 7,3 
0.258  iCl 
C. 235£58 
C.332( 15 
C. 367572 
C. '•06525 
C.'i',3'!,6  7 
C.  '•30^‘i't 
C.51  /AUl 
0.55^358 
0.571216 
0.628273 
C.06523C 
C.  702  If  7 
C.  730U6 
0. 77o iC2 
C .013055 
C.850C16 
C. 686573 

0. 023531 
C. 960688 
C. 997655 
1.0350C2 
1.071759 

1.  108  7 17 
1.155675 
1.  18262  1 
1.219588 
1.256556 
1.293503 
1. 33056C 
1. 367517 
1.505375 
1.551 332 
1.578239 
1.515256 
1.5522C3 
1.589161 
1.626  118 
1.663075 
1. 700C32 
1. 736S89 
1.  7 73957 
1. 8109C5 
1.857661 
1.685818 
1.921 776 


M60  GUN  THERMUCOUPLF  10 
TeUPEKATURE 


21  INCHES  FROM 
ANU  HEAl  FLUX 


BREECH 


PREOICTEO  OY 

HSAI-YIN 

LEE 

IREALTM) 

(NOR.FLUXl 

(REAL.FLUXI 

(NOR.  1 

IREAL.SUF.T ) 

ITri.CO.T) 

tSEC.  1 

(OMLESS) 

IHTU/FTSO  SECl 

1 DPLESSl 

(OFC.F. 1 

(CEG.F.l 

O.OOIOCOO 

0.8395 

1505.5883 

0.1525 

160.8969 

78.80256551 

0.0020000 

0.9208 

1650.3766 

0.2550 

215. 6579 

79. 10233268 

0.O03CC00 

0.8650 

1558. 6551 

0.3169 

259. 5596 

80.75107503 

0. 0C50CC0 

0.7536 

1350. 6657 

0.3509 

267. 8668 

85.2552  J705 

0. 0050000 

0.6257 

1121.5055 

0.3638 

275. 8225 

89, 27332529 

0. C060COO 

0.5907 

893.8086 

0.3618 

273. 7626 

95.23538855 

C.CO  70COO 

0.3818 

685.3901 

0.3599 

267.3031 

101.53891376 

0. OOBCCOO 

0.2796 

501.2017 

0.3316 

257.5553 

107.7C955560 

0. C09CC0C 

0.  1936 

357.0157 

0.3098 

255. 7306 

113-51650288 

C. OICOOOO 

0. 1236 

221.5337 

0.2066 

233. 2537 

118.55552591 

0.  oucoco 

0.0685 

122.6979 

0.2635 

220. 7827 

122.72522558 

0. C12CCC0 

0.0265 

57.5115 

0.2515 

2C8.  881 7 

126.2C00C0C0 

C.C13CCC0 

-0.0051 

-7.5187 

0.22  10 

197.8765 

128.90756935 

0. C15CCC0 

-0.0255 

-55.5355 

0.2026 

187.9522 

1 30.9095  76  2'J 

0.015C00C 

-0.0391 

-70.1121 

0.1863 

179. 1825 

132. 20750150 

0. C16CCC0 

-0.0569 

-85.1269 

0.1722 

171.5612 

133.  1312  1585 

0. Cl 7CC00 

-0.0503 

-90.1656 

0.1600 

165.0278 

133.53167289 

0. 018CC00 

-0.0505 

-90.5995 

0.1598 

159.5080 

133.575033C9 

O.C19COOO 

-0.0585 

-86. 8935 

0.1511 

155.8290 

133. 33976926 

0. J2CCU00 

-0.0551 

-00.8193 

0.1339 

150. 9312 

132.895  7 7680 

■7.021CC00 

-0.0509 

-73.3955 

0.12  73 

1 57.  6772 

•32.29869263 

0.O22CCOO 

-0.0365 

-65.5523 

0.1228 

155.9571 

131.6CCCC0C0 

0.0230000 

-0.0321 

-57.5369 

0.1185 

152. 6723 

130.83766159 

0. 0250000 

-0.C279 

-50.0565 

0.1150 

150.7379 

130.05207950 

0.0250CCO 

-0.0251 

-53.1766 

0.1119 

139. 0626 

129.23622713 

0. U26CCC0 

-0.0206 

-37.0095 

0.1092 

137.6590 

128.53685019 

0.0270CCO 

-0.0176 

-21.5388 

0.1069 

136. 3923 

127.65558013 

0.0280000 

-0.0159 

-26.7007 

0.1058 

135.2788 

126.9001 17C3 

0. 029CCOO 

-0.0125 

-22.3991 

0.1030 

135.2837 

126. 17509350 

0.03COCOO 

-0.0103 

- 18.5266 

0.1013 

133.3900 

125.582950C8 

C.0310C00 

-0.0085 

-15.9803 

0.0998 

132.5860 

125.82560683 

0. C32CCC0 

-0.0065 

- 11.6736 

0.0985 

131. 8635 

125.2COOCOCJ 

0.OJ3OOUO 

-0.0058 

-8.5537 

0.0973 

131.2163 

123. 60858629 

0.0350000 

-0.0031 

-5.5551 

0.0962 

130,6395 

123.05912657 

O.035CCC0 

-0.0015 

-2.7011 

0.0953 

130. 1271 

122.52086579 

C.0360C00 

-0.0000 

-0.0017 

0.0955 

129.6728 

122.02262935 

0.0370CCO 

0.0015 

2. 5000 

0.0937 

129.2632 

121.55335075 

0-0300000 

0.0026 

5.7527 

0.0930 

128.9033 

121.11195267 

0. C39CC00 

0.0037 

6.6538 

0.0925 

128.5662 

120.69729537 

0.05CCC00 

0.0055 

8.1551 

0.0918 

128.2537 

120.30810907 

J. 051CCC0 

0.0051 

9.1705 

0.0912 

127.9213 

119.95292268 

0.C52CC00 

0.C055 

9.6319 

0.0905 

127.5851 

119.60000000 

0.0530000 

0.0053 

9.5859 

0.0899 

127.2172 

119.27729167 

C.C55CCC0 

0.0059 

8. 6978 

0.0391 

126.8066 

118.97250901 

0.05500U0 

0.0050 

7.2565 

0.0882 

126.3397 

118.68262358 

0.0560000 

0.0029 

5.1763 

0.0872 

125. 8062 

118.50589326 

0.05  7CC00 

0.0015 

2.5992 

0.0861 

125.  1985 

118.13591670 

0.058C000 

-0.0005 

-0.7C61 

0.0858 

125. 5119 

117.87221030 

C. 059CCC0 

-0.0075 

-5.3552 

0.0835 

123.7561 

11 7.61020775 

0.05CCCCO 

-0.0056 

-8.3012 

0.0819 

122.9038 

117.35637597 

J.O510COO 

-0.0069 

-12.5387 

0.0802 

121.9919 

1 1 7.07733556 

C. C52CCC0 

-0.0093 

-16.6090 

0.07P5 

121.0208 

1 16.8COQOOCO 
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I.SS3T33 

0.0530000 

-0.0115 

-20.6569 

0 .0765 

120.0095 

116.51169219 

1.'59‘j6S0 

0. 059CCC0 

-0.0136 

-29.9296 

0.0795 

1 18.9596 

116.21026603 

i.032E97 

0.0550CCQ 

-0.0155 

-27. 7699 

0.0726 

1 17. 9099 

115.39920669 

1. C6S6C5 

0. C56CC00 

-0.0170 

-30.590H 

0.07C6 

1 16.8606 

115.56270890 

106  it  2 

3. 057CCC0 

-0.0182 

-32.6399 

0.0688 

1 15. 8975 

1 15.21572882 

i.  1A3‘  19 

0.058CCC0 

-0.0190 

-33.9759 

0.0670 

1 19. 8858 

1 19.85900933 

2.  ISC-iTt 

Q. 059CC0J 

-0.C192 

-39.9966 

0.0653 

1 13.9991 

119.979091C9 

2.217933 

0.06C0CC0 

-0.0191 

-39. 1E82 

0.0638 

113. 1885 

119. 09306290 

2.259391 

C. 061CCC0 

-0.0185 

-33.0787 

0.0625 

112.9817 

1 1 3.698926  15 

2. 291 39H 

0. 062CCOO 

-0.0179 

-31.2390 

0.0619 

111.8818 

113.3CC3COCO 

2.328305 

0.O63C00O 

-0.0161 

-28.7829 

0.0605 

111.3918 

112.90001591 

2.365262 

0. 069CCCO 

-0.0199 

-25.8693 

0.0598 

111.0099 

1 12.50238  769 

2.902220 

0. 065CC00 

-0.0126 

-22.6726 

0.0593 

1 10. 7258 

112.11252006 

2.939177 

0. 066CC0O 

-0.0108 

-19.9252 

0.0539 

1 10.  5267 

1 1 1.73259289 

2.976139 

0.O67CC00 

—7.0091 

-16.3580 

0.0596 

1 10. 3920 

111.36635320 

2.513091 

0.0600COO 

-0.0077 

-13.7137 

0.0085 

110.2968 

111.01690959 

2.55 JC^e 

0.C69CCC0 

-0.0065 

-11.7275 

0.0583 

1 10.2122 

110.68959671 

2. 587CC6 

0.07C0C00 

-0.0C59 

- 10.6115 

0.0581 

110. 1071 

110.37157113 

2.623S63 

0. 0710000 

-0.0059 

-10.5372 

0.0579 

109.9999 

110.077199  19 

2.660920 

0. 072CCC0 

-0.0065 

-11.6178 

0.0579 

109.7091 

109.30030000 

2.69  76  77 

0.0730000 

- 7.0077 

-13.8892 

0.0567 

1C9. 3599 

1C9. 53790973 

2. 739835 

0.C79CCC0 

-0.0096 

-17.2932 

0.0558 

108.8823 

109.28579062 

i. 711152 

0.075COCO 

-0.0121 

-21.6612 

0.0597 

108.2665 

109.09069932 

2.306  799 

0. 076CC00 

-0.0199 

-26.7011 

0.0533 

107. 5151 

108. 79689356 

2 . 895  706 

0.077CCOO 

-0.0178 

-31.9890 

0.0517 

106.6961 

108.59895963 

2 .882663 

0. C78CC00 

-0.0206 

-36.9667 

0.0999 

105. 6999 

108.29171788 

2. 9196 2 1 

C, C79CCC0 

-0.0223 

-90.997C 

0.0981 

109. 7166 

108.02053902 

2.956578 

0. OdCCCOO 

-0.0291 

-93. 1290 

0.0969 

103. 7867 

107.73207721 

2.993535 

0. 081COCO 

-0.0238 

-92.6269 

0.0999 

103.0010 

1C7. 92987919 

3 .030992 

0. 082CCC0 

-0.0215 

-33.5100 

0.0939 

102.973  7 

107. ICCOOOCO 

3.067950 

0. C83CCOO 

-0.0167 

-29.0692 

0.0937 

102.3392 

106.76167727 

3- 1099C7 

0- C89CCC0 

-0.0C89 

-15.8997 

0.0999 

102. 7200 

106.91733289 

3.  191369 

0. 085CCC0 

0.0022 

9.0196 

0.0963 

103.  7685 

106.08051697 

3.173321 

0.0860000 

0.0168 

30. 1 193 

0.0957 

105.6030 

105. 76611295 

3.215278 

C.087CC00 

0.0397 

62. 1950 

0.0598 

1C8.3163 

105.99529511 

3. 252236 

0.0880000 

0.0553 

99.1291 

0.0615 

1 11.9972 

105.29228652 

3. 289193 

C. C89CC00 

0.0776 

139.0999 

0.0699 

116.9528 

105.18932393 

3.  326  15C 

0.09CCCCO 

0.0998 

178. 8266 

0.0796 

121.6722 

105.  19t'9509S 

J.  363  1C7 

0.0910000 

0. 1191 

213.95S9 

0.0900 

127. 2828 

105.36921092 

3.900C65 

0.09200CQ 

0.1317 

236.0913 

0.1001 

132.7977 

105. 7COOCOC1 

2.937022 

O.C93CCCO 

0.1329 

237.2896 

0.1085 

137.2516 

106.21821278 

3.  9739  79 

C. C99CCC0 

0.1191 

209.5369 

0.1129 

139.6259 

106.9  1338298 

2.510936 

0.C950C0C 

0.0678 

121.5869 

0.1109 

138.2573 

107.  757539  78 

3.597693 

0.0960000 

-0.0180 

-32.2719 

0.0969 

130.9886 

1C8.689980C8 

3.589651 

0.097CCC0 

-0.1579 

-283.0999 

0.0672 

1 19.9929 

109.58959982 

3.6218C8 

0. C98CCC0 

-0.  3699 

-662.SB25 

0.0195 

86.6393 

110.28138952 

3.658765 

0.099CCC0 

-0.6760 

-1211.6918 

-0.0696 

91.3115 

110.52073128 

3-695722 

0. ICCCCOO 

-1.1030 

-1977.0890 

-0.1959 

-26. 7312 

109.99096599 

TEMPERATURE  AND  HEAT  FLUX 

PREDICTED  BY 

HSAI-YIN 

LEE 

1 NOR  T M ) 

( REALTMI 

(NUR.FLUX) 

1 RE AL. FLUX) 

(NOR.  ) 

(REAL. SUE. T) 

CTH.CO.T) 

(OMLESSI 

(SEC.) 

(OMLESS) 

IRTU/FTSO  SEC) 

IDMLESS) 

(DEG.E. ) 

(CEG.E.) 

M60  GUN  THERMOCOUPLE  10  21  INCHES  FROM  BREECH 
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M60  GUN  ThERMCCOuPLE  7 15. U INCHES  FROM  BREECH 
(U)RE  SURFACE  IfMPFRATURE  *ND  MEAT  FtUX  PROGRAP 
NUPOEF  UF  BID  COFF.  TO  HE  FITTEO  >-  10 


TIME 


TEMPERATURE 


O.OIOOOOOCOO 

O.02OOOC0000 

O.O3OOOCO00O 

O.C^OOOCOOOO 

0.C500000000 

G.C600000000 

O.C700000000 

0.080000COOO 

O.C9000UOOOO 

O.ICCOOOOCOO 


TIME 

BORE 


CF  DATA 

SHIFTED 

3YISEC)  -0 

FAOlUS 

IHT.)  = 

0.01250 

R RADIUS 

(HT.)* 

0.0AA17 

TC  THERMOCOUPLE 

DISTANCE 

lie. 6000000000 
12A.5COOOCOOOO 
123.10OO0C0OOO 
119.6000000000 
1 16.9C000C0C00 
IIA.IOOOOCOOOU 
1 12.5CCC00U000 
11O.3COOOO0OOO 
107.60COOOOOOO 
1C6.BCOOOCOCOO 


EORE  TC  THERMOCOUPLE  DISTANCE  (FT»=  0.001830 

INIIIAL  thermocouple  TEMPERATURE  IF.)=  7R.8000 

INITIAL  GAS  TEMPERATURE  (F.)*  <r.A937 

THERMAL  niFFUSIVITY  (FTSO/SEC)  = C.0OCIO3O7 

THERMAL  CCNDCCT  I VITYI 8TU/FT.SEC.F . 1=  C. 00555555 

NUMBEF  CF  TIME  TEMPERATURE  PAIRS  (SEC. ,F. 1 = 10 
NUMBEF  UF  fid)  CCFF.  TO  BE  FITTED  = 10 


TIME 


TEMPERATURE 


0.0120000CCO 

O.O223OCO00O 

O.C320COOOOO 

O.CA2OOOOO0O 

0.0520CCCCC0 

0.O62OCC0O0O 

0.0720000000 

O.C82COOCOOO 

O.C9200COOOO 

0.1C20000000 


1 10.6COOOOOOOO 
12A.50COOCCOOO 
123. lOOOOCOUOO 
1 19.6CCOOCOOOO 
1 16.90UC0C00O0 
1 lA. lOCCOOOOOO 
1 12.5COOOCOCOO 
1 10.3000000000 
107.6CCCOCOOOO 
IC6.80O00C0OO0 


IlME  CF  DATA 

SH 

IF  TED 

BY! SEC)  -0.002000 

COEEF  ICIEMS 

UF 

Bd  ) = 

0.A2998293D 

01 

I 

1 

COEFFICIENTS 

OF 

8111  = 

-0.237957C30 

02 

s 

2 

COEFF  ICIEMS 

OF 

B(I)  = 

0.6A87A208D 

02 

I = 

3 

COEFF  ICIEMS 

OF 

ad  )= 

-0. 10652A20D 

03 

9 

4 

CUEf  F ICIEMS 

OF 

B(  I)  = 

0.112355A30 

03 

I = 

5 

COEFF  ICIEMS 

OF 

Bl  I)  = 

-J. 7759011AD 

02 

1 * 

6 

COEFF  ICIEMS 

OF 

BI  I)  = 

0.3A7761130 

02 

1 S 

7 

COEFF  ICIEMS 

OH 

Bld- 

-0.971723300 

01 

I * 

8 

COEFF  ICIEMS 

OF 

blll- 

0. 153A3C23D 

01 

I = 

9 

COEFF  ICIEMS 

OF 

Bd)  = 

-0. 10A3A2130 

00 

1 « 

10 
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M60  CUN  rHeuMccuuPL e i lo.o  inchis  frum  breech 

TEMPI  RATURE  AND  HEAT  FLUX 
PHEDICTED  BY  HSAl-YIN  LEE 


{ NURTM ) 

( RE  AL  TM) 

(NL3R.FLUX) 

(REAE-ELUX) 

(NOR.  ) 

(REAL-SUF. T ) 

( TH 

(OMIE  SSI 

IStC.l 

1 UMLESS) 

(BTU/ETSC  SECJ 

( CMLCSS) 

1 OEG.F. ) 

(DEC 

(.030 J77 

O.OOIOOOO 

0.6763 

1106.1522 

0.1116 

138.9277 

78 

0.061555 

C. 6020000 

0.7553 

1235. 3923 

0.1832 

180.2072 

78 

0.092332 

U. C03C060 

0. 7260 

1187.5239 

0.2382 

207.  1 31  3 

79 

C.  123 1C9 

0. 0090000 

0.6538 

1069.9669 

0.2682 

223. 3065 

31 

C.  153E86 

0.  VJ05OCO0 

0. 56-6  5 

926.6396 

0.2336 

231.5973 

83 

C . 189669 

0. 006C600 

0.9779 

781.6665 

0.2885 

239.252 1 

67 

0.  2159M 

0. LC7COCO 

0.3950 

696.0729 

0.286? 

233.0109 

9 1 

0. 296210 

O.OC8CCOO 

0.3211 

525.3098 

0.2791 

229,  1991 

95 

C.  2769':5 

0. 0O9CCOO 

0.2576 

921.3085 

0.2691 

223.8069 

99 

0.307773 

0. OICCCCO 

0.2092 

339.0C59 

0.2575 

217. 5565 

103 

0. 336550 

0.0110000 

0. 1603 

262.1869 

0.2953 

210.9596 

107 

C. 369327 

C. C12CCC0 

0.1298 

209.0755 

0.2330 

209. 3639 

1 10 

0.900  10  5 

0. Cl 3CCG0 

0.0969 

1 57.6791 

0.2212 

157.9933 

113 

C.930E82 

J. 0190000 

0.0790 

121.0127 

C.2101 

191.9790 

115. 

0. 961659 

C.L'lbOCOO 

0.0569 

92.2309 

0.1997 

166.3861 

117 

C. 992936 

0. C 16CCC0 

0.0926 

69.7006 

0.1901 

181. 2333 

119. 

C.  5232  19 

0. 01 70(00 

0.0318 

52.0336 

0.1313 

176. 511  1 

121 

C. 553991 

0. C18CCC0 

O.0233 

38.0919 

0.1733 

172.  1893 

122 

(.589768 

0. 019C000 

0.0165 

26.9727 

0 .1660 

168.2326 

123 

(.615595 

0. C2CCCC0 

0.0110 

17.9899 

0.1592 

169.6001 

12  3 

0.696 ’23 

0. 021CC00 

0.0065 

10.6397 

0.1530 

161. 2599 

129 

C.6771C0 

0. 022CCOO 

0.CO28 

9. 5700 

0.1973 

158.  1691 

129 

C.  707E  77 

O.023CC0O 

-0.  0003 

-0.9981 

0. 1920 

155.3006 

129 

C. 738659 

0.0290000 

-0.0028 

-9.5583 

0.1371 

1 52. 6936 

129. 

C. 769932 

0. 0250CCO 

-0.C098 

-7.8380 

0.1325 

150. 1 780 

129. 

C. 800 209 

0. C26CC00 

-0.0063 

-10.3301 

0.1232 

197. 8935 

129 

(.  830S  86 

0.  02  70C00 

-0.0079 

-12.0536 

0. 1293 

195. 7836 

129, 

0. 861 769 

0. C28CCC0 

-0.0080 

-13.0279 

0. 12C7 

193.8990 

129. 

(.992591 

0. 0290000 

-0.0081 

-13.2779 

0.1179 

192.0720 

129, 

0.9233  18 

G. 03CCCCO 

-0.0079 

-12.8978 

0.1199 

190. 9696 

123, 

C-9590S5 

J.C310CC0 

-0.0072 

-11.8000 

0.1113 

139. 0180 

123. 

0.989E  73 

0.0320000 

-0.0062 

-10.2193 

0.1099 

137. 7265 

123. 

1.01565C 

0,0330000 

-0.CC50 

-8.2111 

0. 1072 

136.5822 

122. 

1.096927 

0.0390000 

-3.0036 

-5.8979 

0.  1059 

135. 5795 

122, 

1.077209 

0. C35CCC0 

-0.0021 

-3.9195 

0.1037 

139.6903 

122. 

1.107982 

0.036CC00 

-0.0006 

-0.9C20 

0.1023 

133.9191 

121. 

1. 138759 

0.037CC00 

0.0009 

1.9986 

0.1010 

133.2281 

121. 

1.  1695  36 

0.038CCC0 

0.0022 

3.6537 

0.0999 

132.6132 

120. 

1. 200319 

0.C390CC0 

0.0033 

5.9931 

0.0983 

132.0999 

120, 

1.231C91 

C.C9CCCtO 

0.0091 

6.7669 

0.0978 

131.5166 

120. 

1.261668 

0. 0910000 

0.0096 

7.5981 

0.0969 

130.9956 

119. 

1.292695 

0.0920000 

0.0097 

7.7916 

0.0959 

130.9685 

119, 

1. 323923 

0.U93CC00 

0.0095 

7.3318 

0.0999 

129.9200 

119. 

1. 359200 

0.099CLO0 

0.0039 

6.3368 

0.0938 

129. 3377 

119. 

1.389977 

C.095CCC0 

0.0029 

9.8079 

0.0926 

128.7128 

118. 

1.915759 

0. 096CCC0 

0.0017 

2. 8256 

0.0919 

128.0907 

118. 

1.996532 

0.097CC00 

0.0003 

0.5012 

0.0901 

127.32C8 

118. 

1.977  309 

0.0980000 

-0.0012 

-2.0329 

0.0886 

126.5572 

117. 

1.503066 

O.C99COOO 

-0.0020 

-9.6292 

0.0872 

125. 7579 

117. 

1. 538E69 

C. 05CCC00 

-0.0093 

-7.1125 

0.0856 

129.9352 

117. 

1.567691 

0.C51CCC0 

-0.0057 

-9. 3323 

0.0891 

129. 1095 

117. 

1. 6009  18 

0.0520000 

-0.0068 

- 11.1238 

0.0G26 

123.2890 

116. 

XO.TI 
,F  . ) 

80037178 

0925<.6C5 

t)68'7979l 
2619T2C8 
9832 1 7N7 
4932A2^<- 
.981313C3 
661 J6957 
30C  73^68 
73930391 
359  37369 
,600 JCOCO 
,99519669 
39130675 
95651390 
6f>96  36  17 
06507221 
17919891 
09762728 
70912166 
1 7920092 
5CC JCCCO 
6901H933 
77016617 
75738966 
66677881 
51112060 
30195876 
09793335 
75756275 
93996213 
ICC'OCCCO 
79539896 
38128553 
0127C981 
69912808 
2 7937668 
92163655 
57390112 
23695929 
91136965 
6CCCCCC0 
300569  72 
01266071 
739871 1 1 
9 6 5 3 69  56 
20201698 
992592  76 
68969069 
9261 3C98 
165  )S9E7 
9COOCOCO 


, i 


1.631155 

0.053CCC0 

-0.C075 

-12.3507 

0.0811 

122.5936 

1 16.629809CC 

1-661573 

0. 055COOO 

-0.0079 

- 12.858  7 

0.0797 

121. 7553 

1 1 6.  355J6521 

1.69<?75C 

C. 055CCC0 

-0.0077 

- 12.5835 

0.  i.'7  85 

121.0865 

1 16.0  729  50  7 ' 

1.  7?  » / 

0. 0560 COO 

-0.0070 

- 1 1.5581 

0. 0/75 

120. 5091 

1 15.  7 8 7 ; 12(  . 

1.  305 

0. 0570C00 

-0.0058 

-9.5691 

0.0/65 

120. 0385 

115.59865560 

1. 7t  ^00/ 

0.0  5 80  COO 

-0.0051 

-6.651  3 

0.0/59 

1 19. 6867 

1 15-  20  /092 15 

l-3l5£5<i 

J.059COGO 

-0.0019 

-3. 091  1 

0.0755 

I19.5tl2 

115.92127132 

1.^56636 

C.06CCC00 

0-0007 

1.0727 

0.0753 

119. 3632 

115.6 18251 39 

1 - E 7 7 5 I 3 

C- 0610CC0 

0.0035 

5.6503 

0.0753 

1 19.38/1 

115.363  35  798 

1-908151 

0.0620000 

0.0065 

10.5055 

0.0756 

1 19. 52J0 

115.1 CCCLCOO 

1-938568 

0- 063CC0J 

0.C092 

15.0566 

0.0760 

119.7515 

113.35157576 

1 .-769  75  5 

J.065CCCO 

o.oiin 

19.2967 

0.0765 

120. 0235 

113.62075965 

;.OOJ‘.23 

0. C65CC0G 

0.0139 

22. 7950 

0.0771 

120. 3310 

113.5 1C251C6 

;- 031 300 

0.  C66CCOO 

0.0155 

25.2162 

0.0776 

120.  6229 

113.22159622 

J-062C77 

0.  J67C000 

0.0160 

26.2372 

0.0730 

120.8533 

113.05558985 

i.C92?‘;5 

0. 0680000 

0.0156 

25.5661 

0.0783 

120.9735 

1 12.9  1 155  729 

2. 123632 

0. C65CC00 

0.0  150 

22.9636 

0.0782 

120. 9350 

112.73770657 

2.  1555C9 

O.C  7CCC00 

0.0112 

18.2652 

0.0777 

120.6675 

112.68113169 

2.185166 

0-0710000 

0.C07C 

11.5050 

0.0768 

120.1913 

112.58718931 

2.215563 

0. 0720000 

0.0015 

2.5333 

0.0755 

119.5133 

112.  5CCCCCCC 

256  75  1 

C-C73CLCO 

-0.0052 

-8.5525 

0.0735 

1 18.3307 

112.51287519 

2.2773  18 

O.Q  75CC00 

-0.0120 

-20. 0992 

0.0700 

1 16.9386 

112.31655053 

2. 308255 

0.075CCCC 

-0.0210 

-35.3797 

0.0677 

115.2512 

112.20357218 

2. 339C73 

0.0  76CCC0 

-0.0295 

-58.1059 

0.0650 

113. 3062 

112. J6531290 

2.  369  850 

0.077CC00 

-0.0376 

-61.5230 

0.0601 

111.1669 

111.39050075 

2 . 50  3 t 2 7 

0. 078CC00 

-0.0556 

-73.0337 

0.0559 

108.925  1 

1 1 1.673  7819b 

2.531505 

0.0  79CCC0 

-0.0500 

-81.0020 

0.0518 

1C6. 701 7 

11  U508902C7 

2.562 102 

0- J80C000 

-0.0528 

-86.5052 

0.0580 

105. 6565 

111.09135355 

2. 552559 

0-0010000 

-0.0522 

-85.5602 

0.0558 

102.935? 

1 10.  72065125 

2 . 523  7 5^ 

0.0820000 

-0.0575 

-77.6121 

0.0526 

101. 7655 

U0.3CCOOOOO 

2.555513 

0. 083CCC0 

-0.0377 

-61.6278 

0.0518 

LOI.3582 

109.83720257 

2.505291 

J.0850CC0 

-0.0223 

-36.  5355 

0.0529 

101.8957 

10'’.  355251  C7 

2-616C68 

0.0050000 

-0.0011 

-1.7851 

0 .U560 

103. 6061 

ioa.r.52  72839 

2 . 656  65  5 

0.0860000 

0.0260 

52.5338 

0.0517 

106.6518 

1C8. 35392925 

2.67762  3 

0.08  7C0C0 

0.0585 

95.5527 

C.0600 

111. 101 7 

1C  7.90859996 

2. 7C85CC 

0. C08CCCO 

0.0956 

155.81 76 

0.0709 

116.9867 

107.651 17322 

2.739177 

0.0850000 

0.  1326 

216.9337 

0.0852 

125.  1 565 

107.28935705 

2.769555 

0.0900000 

0.  1683 

2 76.1552 

0.U993 

132.2773 

IC7. 19191523 

2.80)732 

0.091CCOO 

0.  1983 

325.3562 

0 . 1 1 50 

150. 7586 

107.28555352 

2. 331 5C9 

0. 0920000 

0.2152 

350.5135 

0.1297 

158.6775 

1C7.6CCJ0000 

2. 062236 

0.093CC0O 

0.2076 

339.5535 

0.1508 

155.6895 

1C8. 15315365 

2.893C63 

0. 0950COO 

0. 1666 

272. 5553 

0.1550 

156.9235 

108.95256036 

2.923051 

0.0950C00 

0.0765 

125.9505 

0.1375 

1 52. 8606 

109.9359  3123 

2.955610 

0.C9600CO 

-0.0819 

-135.C010 

0.1121 

139. 1916 

111. 06C0  1889 

2.505355 

0.097CCO0 

-0.3315 

-552. 3027 

0.0610 

1 11. 6550 

112. 18503325 

2.016  172 

0.0980C0O 

-0.7C09 

-1156.5363 

-0.0259 

65. 8502 

113.107581C5 

2.056550 

0. 0990000 

-1.2257 

-2003. 1892 

-0.1611 

-7.9752 

113.52825998 

2.C77727 

0.  lOCOOOO 

-1.9552 

-3100.1900 

-0.3600 

-115.1663 

113.028236C2 

TEMPERATURE  AND  HEAT  FLUX 

PREDICTED  BY 

HSA  I-YI N 

LEE 

INOR  TM) 

( KEALTM) 

(NOR.FI.UX) 

( RE AL. FLUX) 

(NUR.  ) 

(REAL.SUF. T) 

(TH.CO.T) 

lOMLFSS) 

(SEC.) 

(DMLESS) 

(PTU/FTSO  SFC) 

(DMLtSS) 

(DEG.F. ) 

(DEG.F. ) 

M60  GUN  THERMUCOUPLE  7 15.0  INCHES  FROM  BREECH 
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M60  GUN  THeRMOCOUPLE  ^ <3.0  INCHES  FROM  EREECH 
bURE  surface  TEMPERATURE  AND  HEAT  FLUX  PROGRAM 
number  of  'U  I ) CUFF.  TU  PE  FITTED  * 10 


TIME 


O.CIOOOOOOOO 
C.C2C0CC0C00 
0.C3OO00CC00 
O.OAOOOOOOOO 
C.CbOOOOOOUO 
O.C6COOOOOOO 
O.C700000000 

o.ooooccoooo 

0.0900000000 
0.  ICCOOCOCOO 
TIME  OF  DATA  SHIFTED  BYISEC)  -0.002000 
FORE  RACIUS  (FT.)  = 0.01250 

CUTER  RACIUS  (FT.J=  O.C5000 

PURE  TU  thermocouple  DISTANCE  (FT)=  0.001030 

INITIAL  THERMOCOUPLE  TEMPERATURE  (F.)=  78.8000 

INITIAL  OAS  TEMPERATURE  (F.)=  A.A937 

thermal  OIFFUSIVITY  (FTSO/SEC)  = C.CCC10307 

thermal  conduct  IVI TYl BTU/FT .sec .F . )=  0.00055555 

NUMBER  CE  TIME  TEiMPtRATURE  PAIRS  (SEC.tF.)=10 
AUMHEP  OF  Pin  COFE.  TD  BE.  FITTED  = 10 


TIME 


O.O12OOCO00O 
C. 0220000000 
0.032UOOOQOO 
O.OA20000COO 
0.0520000COO 
O.C620000JOO 
0.0  72JOCOOOO 
0.0320000000 
0.0920000000 
0. 1C2O0CCC00 


T I ME 

CF  DATA 

SHIFTED 

BYISEC)  -0.002000 

CUEFE  ICIEMS 

OF 

B ( I ) = 

0.97265  I860 

01 

= 

1 

CUEFF 

1C lENTS 

OF 

BI  n = 

-0.56818786D 

02 

I = 

2 

CCPFF 

IC lENTS 

OF 

H ( I ) = 

0. 150730510 

03 

3 

CUEFE 

lUIENTS 

OF 

0 t 1 ) = 

-0.228639730 

03 

I = 

4 

C TEEF 

lU lENTS 

OF 

3 ( I ) = 

0.216195640 

03 

I = 

5 

C JtFF 

ICIEMS 

OF 

B I I 1 = 

-0. 131914710 

03 

I = 

6 

CCEFh 

ICIEMS 

OF 

H ( I ) = 

0. 52022  1490 

02 

1 = 

7 

(JEFF 

ICIEMS 

OF 

U 1 I ) = 

-0. 128146550 

02 

f = 

8 

(UlFF 

IC I ENTS 

OF 

B ( I ) = 

0.179336800 

01 

I s 

9 

CUEFf 

ICIENTS 

OF 

B ( I ) = 

-0. 108875930 

00 

1 = 

10 

TEMPERATURE 


1A3.2COCOCCOOO 
157.3COOCCOOOO 
151.0000000000 
lAA .4COCOCOCOO 
137.50COOCOOOO 
132.00C00C0000 
1 29.0COOOCOOOO 
125.9000OC0O0O 
122.80O00C0OOO 
1 19.6000000000 


TEMPERATURE 


1A3.2CCOCCCOOO 
157.3COCOCOCOO 
1 51 .CCOOOCOCOO 
1 AA.ACCCCCOOOO 
137.5CCOOCOOOO 
132.8000000000 
129.COUOOCOOOO 
125.9000000000 
122.8COOOCOCOO 
I 19.60CCOCOOOO 


MhO  GUN  THt  RMUCOUPLf  ^ 9,0  INCHl-S  FROM  EREECH 
lEMPEHATURf  AND  HI  AT  ELUX 
PREDICTED  OY  HSAl-YIN  LEE 


1 NOR IM) 

(RE  A1  IM) 

(NUR.FLUX) 

(REAL. FLUX) 

(NOR.  ) 

(REAL. SUE. T » 

(TH 

(UNLESS) 

(SEC.  ) 

(OMLCSS) 

(HTU/ETSy  SEC) 

(OMLESS) 

(OEG.E. 7 

( CEG 

C.03h  777 

0. OCICCCO 

1.5068 

2969.6179 

0.2997 

213.3551 

78 

0.061 555 

0.0020000 

1.6978 

2695. 2939 

0.9155 

302. 6683 

79 

C.G92  3 32 

J. CO  30000 

1.5392 

2517.6782 

0.5171 

357.9219 

80 

0.  123  109 

0. C09CCC0 

1.3339 

2181.0196 

0.5700 

386. 3269 

89 

C.  153E86 

0. 0050000 

1.0973 

1795.6995 

0.5895 

396.9195 

90 

C.  1896  69 

C. CC6CCCU 

0.3695 

1919.0289 

0.5637 

393.2918 

9 7 

C. 21599 1 

J. J07CCC0 

0.6512 

1065. 1228 

0.5616 

381. 3690 

1C6 

C .296218 

0. 008CCOO 

0,9661 

762.9626 

0.5299 

369. 0339 

119 

0.276995 

C.CC9CC0G 

0.3123 

510.7522 

0.9920 

393. 8737 

123 

0. 307733 

C. CICCCOO 

0.1892 

309.9211 

0.9528 

322. 7799 

13G 

C- 338550 

C.CllCCCO 

0.0996 

159. 80 16 

0.9199 

302.0751 

13  7 

C.  36932  7 

0. 012CCC0 

0.0259 

91.5359 

0.3789 

282.6736 

193 

C . 900  1C  5 

0. 0130000 

-0.0223 

-36.9968 

0.3958 

265. 1212 

197 

0.930682 

C. C19CCC0 

-0.0523 

-85.5691 

0.3172 

299. 701 0 

151 

J.  96  16  59 

0. J15CCCC 

-0.0683 

-111.7083 

0.2927 

236. 9999 

159 

C.  9929  36 

0. C160CCO 

-Q.C736 

-120.9982 

0 .2722 

225.9530 

156 

C . 5232  19 

0. Cl  7C000 

-J.0713 

-1 16.6880 

0.2559 

216.9023 

157 

C.  55  1991 

J.OloOCCO 

-0.0690 

-109. 6290 

0.2919 

209. 1290 

157 

C. 589708 

0. 019CCCO 

-0.0536 

-87.7361 

0.2312 

203.3609 

158 

C-615595 

0.02CCCC0 

-0.C921 

-68. 8C98 

0.2228 

198. 8995 

158 

C. 696323 

0.021CCCO 

-0.0306 

-99.9833 

0.2162 

195.3119 

157 

J.o77  ICO 

0.  0220  COO 

-0.0201 

-32.8107 

0.2111 

192.5163 

157 

C. 7C7E77 

0.O23CCCO 

-0.0112 

-16.3339 

0.2068 

190. 2392 

156 

C. 738659 

0. 029C000 

-0.0099 

-7.1563 

0.2032 

188.291 3 

156 

C. 769932 

0. 025CCC0 

0.0003 

0.9769 

0.1999 

186. 5168 

155 

C.8J02L9 

0. C26CCC0 

0.0028 

9.6C83 

C.1967 

189.7935 

159 

C-830986 

J.  02  7CCC0 

0.0039 

5- 5020 

0. 1935 

183.0319 

159, 

0.661 769 

0. 0280000 

0.0022 

3.5769 

0.1900 

181. 1703 

153 

0.892591 

0. C29C000 

-0.0009 

-0.6965 

0.1863 

179. 1 763 

152 

0.923318 

0. 03CQC00 

-0.0090 

-6.5995 

0. 1823 

177.0339 

152, 

C.959C95 

0. 031CCC0 

-0.0089 

-13.68C8 

0.1781 

179.7637 

151, 

0.989873 

0. 032CCC0 

-0.0130 

-21.3395 

0. 1737 

172.3793 

151, 

1.01565C 

0.0330000 

-0.0178 

-29.0511 

0.1691 

169.9022 

150. 

1. 096927 

0. C39CCC0 

-0.0222 

-36.3615 

0 . 1694 

167.3862 

199, 

1.C772C9 

0. 0350C0C 

-0.0262 

-92. 8952 

0.1597 

169. 8685 

199, 

1 . 1079C2 

0. 0360000 

-0. 0296 

-98.3621 

0.1551 

162. 391 3 

198, 

I.  138  759 

0.  03  7CCC0 

-0.0321 

-52.56C5 

0. 1507 

159.9997 

197, 

1.  U9536 

J. 03800CJ 

-0.0339 

-55.375 1 

0.1965 

157. 7199 

197, 

1.200319 

C. C39CCC0 

-0.0397 

-56.7716 

0. 1925 

155.5802 

196. 

1.231091 

0. C9CCC00 

-0.0397 

-56.7889 

0.1339 

153.6153 

195 

1.261 668 

0. 0910000 

-0.0339 

-95.5290 

0.1356 

151.8353 

195, 

1. 2926h5 

0. 092CC00 

-0.0325 

-53. 1998 

0.1326 

150.2981 

199, 

1.323923 

J. 093C000 

-0.0305 

-99.8289 

0.1300 

198. 8595 

193, 

1. 3592CC 

0. C99CC00 

-0.0230 

-95.7978 

0.1278 

197.6982 

192, 

J. 389977 

0.095CCCO 

-0.0252 

-91.2895 

0. 1259 

196.6170 

192. 

1.915759 

0. 0960000 

-0.0223 

-36.5215 

0.  1292 

195. 7937 

191. 

1.996532 

0. 097CC0U 

-0.0199 

-31.7327 

0. 1229 

195.0C71 

190. 

1.9773C9 

J. C98CCC0 

-0.0166 

-27.  1228 

0.1217 

199.3839 

190, 

1. 5C8C66 

C. C99CC00 

-0.0190 

-22.8698 

0.1207 

193.8h79 

139. 

1.  538669 

0.05CCCCO 

-0.0117 

- 19. 1187 

0. 1 198 

193. 3738 

138, 

1 . 369691 

0.051CCC0 

-0.0098 

-15.977C 

0.1190 

192.9382 

138. 

1.600918 

0.052CCCO 

-O.C083 

- 13.5129 

0.1103 

192.5183 

137, 

XO.T  I 

f.i 

,8003T969 
008<'62  16 
9217EA52 
2 781  7903 
2971 A983 
83lg95 19 
28706625 
93619C53 
2999C929 
86059319 
5 95 1 56  09 
20QOCOOO 
808895  79 
9201  16C3 
122329  1 1 
02966  789 
26331105 
99627789 
19923329 
11192795 
7332C9'>  1 
3CCOOOCO 
7C737797 
05693060 
382o3J92 
7C75C99  1 
09599980 
903J6C61 
781376C2 
17798532 
50599730 
CCOOOOCO 
9 12551  72 
31695999 
20757930 
58016928 
9320  7959 
262339  79 
57157537 
86187822 
136  58279 
900OC00O 
65712997 
91313763 
17911192 
99979732 
73013289 
03955198 
36153299 
71379553 
09295266 
5CCCC000 


i 


1.631  1^5 

0. 053CC00 

-0.0072 

-11.7523 

0.1175 

192.0995 

136.939RS5B1 

1.661SI3 

C, C59CCOO 

-0.C065 

-10.6851 

0. 1167 

191.6511 

136.39660628 

1.69275U 

■J.055CCC0 

-0.0063 

- 10.2629 

0.1158 

191. 1 761 

135.88396353 

1.  Ul'.  2 7 

0. C56CC00 

-0.0069 

-10.9069 

0.1198 

190.6627 

135.39959839 

l.lbi  TC-i 

C.  05  7CCC0 

-O.CC6.7 

-11.0122 

0.1138 

190.1072 

139.92609099 

1. J8SCB2 

0.  056CCO0 

-U.0073 

-11.9555 

0.1127 

139.5111 

139.97569799 

1.  IJSf 

0. C590CC0 

-0.0080 

-13.1020 

0.1115 

138.8799 

139.09059593 

1.9^6636 

0. C6CCCC0 

-0.0088 

-19.3193 

0.1103 

138. 2209 

133.61788259 

1 . d 7 7 < 1 3 

0, 061CC00 

-0.C095 

-15.9597 

0.1090 

137.5951 

1 33.2051  1968 

1. 9Cb  IS  1 

0. J620C00 

-0.0100 

- 16.9189 

0.1078 

136.866J 

132.8COCCJOO 

l.93B'6H 

0. 063CCCO 

-0.0109 

-17.C9C9 

0.1066 

136. 1969 

132. 90068930 

1 . 969  (St) 

0. C69CCC0 

-0.0106 

- 17.9039 

0.1053 

135.5993 

132.00531616 

; .OOJ523 

0. 06.5CCC0 

-0.0106 

-17.3155 

0.  1092 

139. 9365 

131.61955532 

C3 1 3CC 

0. C66CCC0 

-0.0103 

-16.8187 

0.1031 

139. 3678 

131.2265779C 

;.062C77 

0. C67CC00 

-0.0097 

-15.9929 

0. 1022 

133.8996 

1 30.89209069 

i.Q92£6S 

0. 0600000 

-0.0090 

-19.7517 

0. 10 13 

133. 3850 

130.96152121 

1236  32 

0. C69CCC0 

-0.0082 

-13.3995 

0.1005 

132.9727 

130.08592073 

i.  159SC9 

0.07CCCC0 

-0.0072 

-11.8969 

0.0999 

132. 6071 

129.71635C9? 

2.  18S  IfcC 

0. C 71CCC0 

-0.0069 

-10.9056 

0.0993 

132.2789 

129. 3590C162 

2.215963 

0-C72CCOO 

-0.0056 

-9. 1790 

0.U987 

131.9723 

129.CC00C0C0 

2. 296791 

0.073CC00 

-0.0051 

-8.3299 

0.0981 

131.6739 

128.65526990 

i.ZllilP 

O.C79CCOO 

-0.0099 

-7.9883 

0.  J97o 

131-3610 

128.3203S732 

2 . 30  0 7 9 8 

0.C75CCC0 

-0.0051 

-8.2866 

0. 0969 

131.0159 

! 27.99552235 

3 39  0 7 3 

0. C76CCC0 

-0.005  7 

-9.2999 

0.0962 

130. 6190 

127.6802-^956 

2.36.9  190 

0.  0 7 7CCC0 

-0.0067 

-11.0318 

0.0953 

130.  1593 

127.37360650 

2.90062  7 

0.  07(30000 

-0.0032 

-13.9953 

0.099  3 

129. 6107 

127.07909017 

..931909 

0.0  79COCO 

-0. 0 loo 

- 16.9026 

0.0931 

128.9891 

126.77997573 

2.962  11)2 

0. 08CCCOO 

-0.0120 

- 19.6820 

0.0913 

128.  2796 

126.98792815 

2.  992  959 

0.O81OCOO 

-0.0190 

-22.9719 

0.0709 

127. 5129 

126. 15517657 

2.523736 

0. OH  20  000 

-0.0158 

-25. 8732 

0.0839 

126. 7122 

125. 9C0  )00C0 

2 . 559  5 13 

J.083COOO 

-0.0171 

-27.9132 

0.0875 

125.9101 

125.55996C12 

2.5d5291 

0. 089CCC0 

-0.0175 

-23.5920 

0.0861 

125.1835 

125.291 73369 

2. 616C6H 

J. 085CU0C 

-0.0167 

-27.3796 

0.0850 

129. 5719 

129.9  7596275 

2.696695 

1.0860CCO 

-0.0195 

-23.7975 

0.0892 

129. 1513 

129.65260591 

2.67762  ) 

0.Q87CCC0 

-0.0107 

-17,5217 

0.0839 

123.9929 

129.32375055 

2. 7069C0 

0.038CCC0 

-0.0052 

-8.9386 

0.0892 

129. 1570 

123.99335278 

2.  739  17  7 

C. 089CCCO 

0.0020 

j. 1999 

0.0352 

129.6832 

123.66733266 

2. 769959 

0.09CCC00 

0.0102 

16.6198 

0.0868 

125.5731 

123.35397992 

2.300732 

0. 09  1000(0 

0.0185 

30.3256 

0.0890 

126. 7683 

123.06103571 

2.d3l 5C9 

0.C920000 

0. 0256 

91.8599 

0.0915 

128.1239 

122. 3C00C0C0 

2.062  266 

0.C93CC00 

0.C290 

97.5118 

0.0939 

129. 3706 

122.57977337 

2.893C63 

0. C99CCC0 

0.0257 

91.9786 

0.0952 

130.0785 

122.90733237 

2.923891 

0. 095CCCC 

0.0110 

17-9579 

0.0993 

129.5999 

122.28952910 

2.959618 

0.0960000 

-0.0210 

-39.323C 

0.0095 

127.0103 

122.20996088 

2.905395 

0.  09  70000 

-0.C780 

-127.6983 

0.0789 

121.0395 

122. 1 9667570 

3.016  17? 

0.0980CUC 

-0.1702 

-278.998C 

0.0573 

IC9.9590 

122.07093581 

3-09o<50 

C. C99CC00 

-0.3103 

-5C7.5151 

0.0236 

SI. 5391 

121.90966125 

3.C77727 

0.  ICCCCCO 

-0. 5190 

-890. 8277 

-0.3296 

62.8978 

121.55763221 

TEMPERATURE  AND  FEAT  FLUX 

PREOICTED  BY 

HSA  1-Y(N 

LEE 

(NOR  FMI 

(HEALIM) 

(NUR.FLUX) 

(REAL. FLUX) 

(NOR.  ) 

(REAL.SUF.I  ) 

(TH.CO.T) 

(DMUESS) 

(SEC.) 

(DMLESS) 

(BTU/FTSO  SEC) 

( CMLESS) 

(OEG.F. ) 

( CEG.F. ) 

M60  GUN  THERMOCOUPLE  4 9.0  INCHES  FROM  BREECH 
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M60  GUN  THERMUCOUPLE  10  21  INCHES  FROM  ERFECH 
BORE  SURFACE  TEMPERATURE  AND  HEAT  FLUX  PROGRAM 
NUMBER  OF  B(I)  COFF.  TO  BE  FITTED  » 10 


TIME  TEMPERATURE 


0.0100000000  126.2COOCCOOOO 

0.C2000U0U00  131.60COOCOOOO 

0.03O00C0000  12A.2COOCCCOOO 

O.OAOUOOOOOO  119.6000300000 

O.C500000000  116.8000000000 

O.C6COCOOOOO  113.3000000000 

O.C700000000  1C9.BOOOCOOOOO 

0.0800000000  107.10000C0000 

O.C9COOOOOOO  105.7C00000000 

C.1C00003000  ICA. 2000000000 

TIME  CF  DATA  SHIFTED  BV(SEC)  -0.002000 
BORE  RADIUS  IFT.)  = 0.01250 

CUTER  RADIUS  IFT.)=  0.03562 

EORE  TC  THERMOCOUPLE  DISTANCE  |FT)=  0.001670 

INITIAL  THERMOCOUPLE  TEMPERATURE  (F.)=  78.8000 

INITIAL  GAS  TEMPERATURE  (E.)=  A.A«>37 

IHERMAL  DIEEUSIVITY  (FTSG/SEC)  = 0.C0C10307 

THERMAL  CCNOUCTIVITYI BTU/FI.SEC.E.)=  0.00555555 

NU''BER  CE  TIME  TEMPERATURE  PAIRS  ISEC.tF.)  = 10 
NUMBER  OE  DID  CUEF.  TO  BE  FITTFO  = 10 


O.O12OCUO000  126.2COOOCCOOO 

0.0220000000  1 ?1.6COOOOOOOO 

0.C320000000  12A.2COOCCOOOO 

O.OA200COCOO  1 19.60COOCOOOO 

0.0520000000  1 16.BCOOCCOCOO 

0.0620000000  1 I3.300C0C0000 

0.0720000000  1C9.3COOOOOOOO 

0.08200C0000  1C7. lOCCOCOOOO 

0.C92000J000  105. 7000000000 

O.1O2COCO0OO  10^.2000000000 

TIME  CF  DATA  SHIFTED  BYISECJ  -0.002000 
COEFFICItNTS  OF  BII)=  0.^93585860  01  1=1 

COEFFICIENTS  OF  Btl)=  -0.23331221D  02  1=2 

COEFFICIENTS  OF  B(I)=  0 . 5 3 A635'’.  3D  02  1 = 3 

CJtFFICIFNTS  OF  i3  ( I » = -0.6SB6A89AD  02  1 = A 

COEFFICIENTS  OF  HII)=  0.578930860  02  1=5 

COtEHlCIENTS  OF  S(I)=  -0.31318033D  02  1=6 

COEFFICIENTS  OF  BII)=  0.11016609D  02  1=7 

COEFFICIENTS  0FBI1)=  -0.2A257631D  01  1=8 

COEFFICIENTS  OF  fl(I)=  0.303290150  00  1=9 

COEFFICIENTS  OF  BII)=  -0. 16 A 136 78D-0 1 I =10 


APPENDIX  D THE  CASE  OF  OSCILATORY  SURFACE  TEMPERATURE 


Consider  a slab  with  a sufficient  thickness,  I,  such  that  when  a surface 

is  subjected  to  a periodic  surface  temperature  variation  with  a frequency 

w the  other  surface  is  held  at  the  initial  temperature  T . If  properties 

o 

are  assumed  constant  the  governing  equation  for  the  problem  can  be  written 


3v  3^v 

— = «— 

3t  3x 


with  initial  and  boundary  conditions  as 


v(x,  o)  = 0 

v(H,  t)  = sin  wt 

v(o,  t)  = 0 


where  v = (T  - T )/(T  - T ) and  a is  the  thermal  diffusivity. 

o max  o 

The  solution  of  the  problem  according  to  Carslaw  and  Jaeger  [4]  can 
be  written  as 


V = ^ 2 (-D(-l)V 

i n=l 


2 2 2 

an  IT  X/i.  . , ,, 

e sin  wXdX 


we  have 


sin  bx  dx  = e**(a  sin  bx  - b cos  bx)/(a^  + b^)  (6) 


PART  in 


■RANSTENT  SURFACE  HEAT  FLUX  AND 
.TURE  ON  A HOLLOW  CYLINDER 


INTRODUCTION 


In  the  study  of  transient  heat  transfer  many  efforts  have  been  made 
on  the  so-called  "inverse  problem"  [1,2]  where  a surface  heat  flux  and 
temperature  is  to  be  predicted  by  the  measured  data  at  some  location  in- 
terior to  a body. 

In  the  previous  works  [1-6]  the  solution  is  represented  in  either  an 
Integral  form  after  some  manipulation  of  the  contour  Integral  from  the 
— inverse  transform,  or  in  a series  form  after  the  expansion  of  the  solution 

for  small  and  large  times.  Using  Laplace  transformation  Chen  and  Thomsen 
[6]  Introduced  a polynomial  in  terms  of  the  error  function  to  represent  the 
response  of  thermocouple  measurement  and  the  inversion  is  accomplished 
for  any  transient  surface  heat  flux  at  the  inner  surface  of  a cylindrical 
tube.  However,  their  Inversion  solution  is  valid  only  for  a short  duration 
due  to  the  asymptotic  expansion  of  the  modified  Bessel  function  in  the  in- 
verse Laplace  transform.  In  this  study  an  exact  solution  obtained  from  the 
Inverse  Laplace  transform  by  the  convolution  method  is  given  for  the  case 
of  hollow  cylinder.  The  solution  is  valid  for  both  constant  and  variable 

I 

j heat  flux  and  for  both  short  and  long  time  duration. 

I ANALYSIS 

I 

Consider  a long  hollow  cylinder  with  sufficient  wall  thickness  such 
that  the  outer  surface  temperature  has  a negligible  response  when  the  inner 
surface  is  exposed  to  a thermal  pulse  of  a transient  process.  This  condi- 
tion considerably  simplifies  the  theoretical  analysis  as  the  outer  boundary 
may  be  assumed  to  be  infinite,  and  only  one  interior  probe  of  the  cylinder 
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I 

1 


I 

I 


Is  required  In  the  experimental  measurement.  The  material  of  the  cylinder 
is  considered  to  be  homogeneous  and  isotropic  with  constant  thermal  dlf- 
fuslvlty,  a.  Let  and  be,  respectively,  the  inner  and  outer  surface 
radii.  Rj^  the  radius  of  the  probe  location  and  t the  dimensionless  time. 
If  the  temperature  of  the  cylinder  is  initially  uniform  at  T^,  the  mathe- 
matical problem  governing  the  temperature  T,  may  be  written  as 


9®.  = + i ii 

9t  ° . 2 r 3r 
dr 


1 < r < r * » 
o 


(1) 


0 (r,  o)  = 0 


(2) 


0 (“>,  t)  = 0 


(3) 


1 


1 

( 


0 t)  = f(t)  1 < < “>  (4) 

2 

where  0 = T - T , r = R/R  , t = ot/R.  , and  f(t)  is  the  interior  tempera- 
o 1 i 

ture  response  of  the  thermocouple  measured  at  r = r^  at  the  dimensionless 
time  t.  The  problem  is  to  predict  the  surface  temperature  9 (1,  t)  and 
heat  flux  per  unit  area 

q = - (K/R,)O0/3r)  I (5) 

^ 'r  = 1 

where  K is  the  thermal  conductivity. 

The  problem  can  be  solved  by  Laplace  transformation.  Let  the  trans- 
formation be 


When  9 satisfies  the  Dirichlet's  condition  the  temperature  function  9 is 
recovered  by  inversion  of  the  Laplace  transformation  as 


0(r,  t) 


1 

2iri 


-c+i“ 


9e^^ds 


(7) 


where  c is  a suitable  positive  value.  Equation  (1)  and  (2)  under  trans- 
formation (6)  becomes 


d^9  . 1 ^ 

.2  r dr 
dr 


s9 


which  has  a solution  of  the  form 

9 = AI  (pr)  + BK  (pr)  (8) 

o o 

where  I and  K are  modified  Bessel  functions  of  the  first  and  second  kind 
o o 

1/2 

with  p = (s)  . With  the  boundary  conditions  (3)  and  (4),  (8)  becomes 

9 = f(s)  [K^(pr)/K^(pr^)]  (9) 

where  f(s)  is  the  Laplace  transform  of  the  boundary  condition  (4). 

The  temperature  response  measured  at  r = r^^  can  be  expressed  by  a 
polynomial  or  numerous  other  suitable  functions.  In  the  present  analysis, 
for  reasons  to  be  explained  later,  f(t)  will  be  represented  as 

N 

f(t)  = E b 
n*l 

1 At 

If  we  choose  F^(t)  = y-  e and  Fn(t  - t)  being  any  arbitrary  function 


"A 


(t)  Fn(t  - x)dT 


(10) 
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r 


depending  on  n,  for  example  (t  - t)  etc.  then  the  Laplace  transform  of 
Eq.  (10)  gives 


f(s)  = 2 (pr  )Fn(s) 

no  1 

n=l 


where  Fn(s)  is  the  Laplace  transform  of  Fn(t).  Substituting  Eq.  (11)  into 
Eq.  (19)  we  have 


e(s,  r)  = ^ b K (pr)Fn(s) 

n o 


It  is  noted  that  the  K (pr, ) in  Eq.  (9)  has  been  cancelled  by  this  sub- 

o 1 9 

r..  ^ 


stitution  which  explains  the  choice  of  F 
The  inversion  of  Eq.  (12)  gives 


(t)  = e~^ 


in  Eq.  (10). 


N ^t 


(t,  r)  = iC  ^ ® 

n=l 


At  surface  r ■ 1 


N ^t 

Kt.  D-i:  \ J 

n=l 


The  temperature  gradient  and  hence  heat  flux  at  surface  is 


ii  V s ^ 

r=l  ,^1  " i 


some  examples  of  Fn(t  - t)  function  and  the  representation  of  the  thermo- 
couple response  are: 


' V ' 


■'•V,  VXi,' 


If  Fn(s)  = —2 2 

s + n 


, " ^ r 1 At  1 

then  f(t)  = E \ / 2t  ® n 

n=l  •'o 


— e - sin  n(t  - 


If  Fn(s)  * ~2  Y 
s +n 


then  f(t) 


N -t 

-2:»J  i 

n=l  •'o 


cos  n(t  - T)dT 


If  Fn(s) 


(s  + a) 


then  f(t) 


N 


“ 4t  . .n-1  -a(t-T) 

_e - t)  e 

2t(n  - 1) 


If  Fn(s)  = s 


-(n+1/2) 


A e 2^(t  - T 

■ li  i 2, 


n-l/2 


(2n  - 1)A 


If  Fn(s)  =»  — 
s 


then  f(t) 


N J 

-Zb  f 

n=l  " •'o 


t 


2 (n-1): 


“•  V ■»* 


DISCUSSION 


In  practical  application  one  normally  has  a measurement  of  an  Interior 
thermocouple  response  preferably  near  the  heating  surface.  The  temperature 
measured  by  the  thermocouple  can  then  be  fitted  with  a suitable  polynomial 
given  in  Eq.  (16)  through  (20).  From  Chen  and  Thomsen's  (6)  work  the  mini- 
mum degree  of  pol3momlal  N should  be  about  10  for  a transient  heating  such 

as  in  gun  barrel  and  engine.  The  coefficients  b obtained  from  the  fitting 

n 

is  then  substituted  into  Eq.  (14)  for  the  surface  temperature  and  Eq.  (15) 
for  the  surface  heat  flux. 
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